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1.0  INTRODUCTION 

1.1  General 

The  knowledge  of  the  hydrologic  conditions  which  exist  during  low  flow  conditions  can  be  of 
primary  importance  for  some  watercourses.  For  example,  when  analysing  water  quality  conditions, 
the  low  flow  characteristics  of  a  watercourse  are  a  major  concern  to  both  the  landowner  and  the 
Ontario  Ministry  of  the  Environment.  Other  uses  of  low  flow  information  may  include  the 
following: 

1)  Instream  pollutant  analyses  (point  and  non-point  sources) 

ii)  Reservoir  design  Qow  flow  augmentation) 

iii)  Environmental  appraisals 

iv)  Feasibility  of  small  hydro  developments 

v)  Water  supply  and  evaluation  for  water  taking  permits 

vi)  Base  flow /groundwater  recharge  and /or  contamination  analysis 

vii)  Stream  fisheries  assessments 

viii)  Analyse    effects    of   changes    in    watershed    on    low    flows    (e.g.    deforestation, 
urbanization) 

ix)  Agricultural 

The  identification  of  low  flow  characteristics  within  a  watercourse  is  most  easily  accomplished 
using  continuous  hydrometric  data  recorded  for  the  stream. 

A  primary  source  of  information  describing  drought  conditions  is  the  "Low  Flow  Characteristics" 
maps  which  were  recently  up)dated  by  Cumming  Cockbum  Limited  for  the  Ministry  of  the 
Enviroiunent.  This  update  also  included  individual  reports  on  the  regions  listing  gauging  station 
data,  results  of  non-parametric  testing  of  low  flow  data,  tables  and  graphs  of  low  flow  values  and 
characteristics  for  both  extreme  value  and  flow  duration  analysis.  A  large  data  base  exists  with 
all  of  this  data.  However,  the  effectiveness  of  utilizing  single  station  analyses  is  limited  since  a 
hydrometric  recording  gauge  may  not  be  located  in  the  vicinity  of  the  study  site. 
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Useful  techniques  do  not  presently  exist  for  transferring  this  tnfomnation  to  ungauged  sites.  The 
use  of  appropriate  techniques  (e.g.  station  proration  by  unit  flows,  area,  etc.)  are  limited  by  several 
assumptions,  including: 

•  ignoring  the  effects  of  regulation  or  upstream  storage  flakes,  swamps,  etc.) 

•  assuming  the  watersheds  are  homogeneous  (i.e.  physiographic  characteristics  are  ignored) 

•  assuming  the  climatic  regions  are  homogeneous. 

Cumming  Cockbum  limited  recently  completed  a  preliminary  research  program  which  describes 
the  initial  stages  of  development  of  techniques  to  produce  low  flow  estimates  for  ungauged  sites 
for  the  Southwestern  and  West  Central  regions  in  Ontario. 

The  study  included  a  literature  review  of  similar  relevant  investigations.  In  all,  five  methods  for 
estimating  low  flows  were  examined. 

The  recent  investigations  have  led  to  the  identification  of  several  areas  of  research  to  be  fiirther 
develop)ed  for  estimating  low  flows  for  ungauged  watersheds.  The  present  study  was  undertaken 
to  confirm  the  applicability  of  regional  transfer  to  modify  and  enhance  the  applicability  of  available 
methods,  and  to  develop  further  insight  into  predicting  low  flow  values  at  ungauged  locations. 
The  test  areas  selected  were  the  Central  and  Southeastern  regions  of  the  Province  of  Ontario  (see 
Figure  1.1). 

1.2         Study  Objectives 

The  main  objective  was  to  further  refine  techniques  for  providing  estimates  of  low  flow 
characteristics  for  imgauged  streams  based  on  the  physical  parameters  of  the  watershed  and 
appropriate  meteorological  variables. 

It  is  expected  that  the  technique  could  then  be  further  develojjed  and  adapted  in  order  to  provide 
estimates  of  low  flows  for  imgauged  watersheds  at  other  locations  in  the  Province. 

The  following  points  summarize  the  focus  of  this  investigation: 


B  .       Central      Region 

C         Southeastern    Region 


CBmnimg  Cockburn  Limited 

Consulting  Engineers  and  Planners 


OCl 


STUDY     LOCATION 


F  I6URE     -     1  •! 
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1)  To  test  the  available  methodologies  for  predicting  low  flows  in  the  Central  and  Southeastern 
regions.  To  identify  suitable  techniques  for  application  and  required  research /refinements 
(e.g.  by  comparison  of  regression  and  mapping  techiuques) 

2)  To  develop  an  appropriate  data  base  including  additiorial  parameters  such  as  evaporation, 
groundwater  fluctuations  (well  records),  etc. 
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2.0  LITERATURE  REVIEW 

2.1  General 

This  updated  literature  review  augments  the  previous  search  (Cumming  Cockbum  Limited,  March 
1990)  with  recently  published  information  on  low  flow  characteristics.  A  summary  of  the  relevant 
results  of  the  review  of  the  various  investigations  is  given  in  Table  Zl  and  discussed  in  the 
Summary  section  Z2. 

The  following  sections  summarize  highlights  from  selected  references.  The  reader  should  refer  to 
the  original  reference  before  attempting  to  apply  specific  equations  since  the  units  and  parameter 
definitions  vary  from  report  to  report.  The  intent  here  is  only  to  identify  important  parameters  for 
further  consideration  in  this  investigation. 

Probabilitv  Distribution  of  Low  Flows  (Matalas,  N.C.,  1963) 

This  study  examined  the  fit  of  four  probability  distributions  using  both  the  method  of  maximum 
likelihood  and  moments.  The  Gumbel  and  Pearson  Type  in  distributions  fitted  the  data  equally 
weU  and  were  more  representative  of  the  probability  distribution  than  the  3  Parameter  Log  Normal 
or  Pearson  Type  V  distributions.  Goodness  of  fit  was  based  on  observed  minimum  flow,  the  lower 
limits  of  the  four  probability  distributions  and  the  relation  between  skewness  and  kurtosis  of  the 
low  flow  data.  It  was  recommended  that  future  statistical  studies  of  low  flows  should  lise  the 
method  of  maximum  likelihood  estimates  of  parameters  rather  than  moment  estimates. 

Modelling  of  Low  Flow  Frequencv  Distributions  and  Parameter  Estimation  (Condie,  R.  and  G.  Nix, 
1974) 

The  fit  of  four  theoretical  frequency  distributioris  was  examined  in  this  study.  The  Gumbel  III 
distribution  was  selected  as  the  best  fit  for  extreme  value  analysis  of  low  flows  since  it  produced 
33  acceptable  fits  out  of  a  possible  38  for  the  Canadian  rivers  used  in  this  study.  A  program 
FLOINT  was  written  to  handle  low  flow  values  of  zero  using  a  joint/conditional  probability 
function.  A  discussion  of  the  lower  boundary  limit  being  between  zero  and  the  lowest  recorded 
flow  concluded  that  the  lower  boundary  can  be  less  than  zero  since  it  represents  the  ciirve  of  the 
Line  and  not  the  actual  flow. 
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Hvdroloev  of  7  Day  10  Year  Low  Flows  (Singh,  K.P.  and  J.B.  Stahl,  1974) 

This  paper  discusses  three  pertinent  factors  in  the  broad  headings  of;  1)  basin  factors;  2)  climatic 
factors,  and  3)  sociological  factors.  The  main  discussion  is  based  upon  man's  influence  and 
requirement  for  low  flow.  It  is  suggested  that  basin  input  to  a  7Qio  low  flow  study  includes 
information  on  daily  streamflows,  soil  and  stream  characteristics,  physiographic  and  hydrologic 
diversion  boundaries,  municipal  and  industrial  water  use  and/or  man-made  lakes  and  flow 
regulation.  Again,  it  is  stressed  in  the  conclusions  that  the  magnitude  of  the  drainage  area  is  the 
most  important  factor  affecting  the  natural  7Q,o  of  hydrologically  similar  basins. 

Determining  Streamflows  from  Geomorphic  Parameters  (Osbom,  J.F.,  1974) 

Osbom  examined  the  many  interconnected  relationships  of  geomorphic  parameters  of  basins  to 
flows.  An  example  of  the  basin  characteristics  interrelationship  is  that  LNTH  =  1.4  (DA)°*,  where 
LNTH  is  stream  length  and  DA  is  area.  It  is  noted  that  the  coefficient  and  power  values  vary  for 
different  geographic  areas. 

The  equations  for  estimating  average  annual  flow  (QAA)  took  the  form  QAA  -  C"  (P)'  QA)''  where 
P  is  average  annual  precipitation  and  the  coefficient  and  powers  are  determined  from  regional 
regressions. 

It  is  noted  that  minimum  or  low  flows  tend  to  be  more  sensitive  to  a  larger  number  of  parameters 
than  flood  flows.  Assuming  linear  flow-frequency  relationships,  a  parameter  Q7LIP  was  developed 
to  relate  the  relative  maximum  size  of  low  flow  in  one  stream  to  another.  This  value  was  then 
regressed  with  basin  parameters  to  establish  IQ^. 

.  ^    PAH)''  Q7L1P 
^  3OO(/0 

where  B  is  a  coefficient 

H  is  basin  relief 
P  is  slope  of  recurrence  interval  graph 
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The  P  parameter  also  represents  a  measure  of  the  stability  of  the  natural  basin  storage  from  which 
the  low  flow  is  derived.  The  value  7Q2  was  then  used  as  an  index  flow  to  predict  other  recurrence 
interval  flows  as: 


where  P  is  slope  and  C  is  a  coefficient.  For  example, 

log  7Q20  =  log  (P(7Q2))  -  log  C 

Relationships  were  also  found  between  peak  flows  and  low  flows.  It  was  determined  that  if  basins 
were  of  the  same  average  annual  runoff  then  you  assume  C  is  essentially  the  same  and  examine 
the  peak  flows,  if  they  are  different,  then  low  flows  must  be  different  and  the  basins  are  not 
hydrologically  similar. 

Combining  Estimates  of  Lx)w  Flow  Characteristics  of  Streams  in  Massachusetts  and  Rhode  Island 
(Tasker,  G.D.,  1975) 

This  report  presents  a  statistical  analysis  of  baseflow  measurements  and  a  discussion  of  the 
accuracy  produced  by  varying  the  sample  size.  It  was  concluded  that  empirically  estimated  values 
of  variance  of  the  estimated  7Qjo  low  flow  from  its  true  value  indicate  that  more  than  six  or  eight 
baseflow  measurements  add  little  to  the  confidence  with  which  such  an  estimate  is  made.  As  well, 
the  number  of  baseflow  measurements  needed  to  meet  a  spsedfic  accuracy  of  low  flow 
characteristics  can  be  sigjuficantly  reduced  by  considering  prior  knowledge  of  the  low  flow 
characteristics  quantitatively.  The  formula  is: 

^        1  +  {Sr^lSb^ 


where  q     is  logarithm  of  estimated  low  flow  characteristics 

q,    is  log  Qr  flow  estimated  by  regression  of  basin  parameters 
qj,    is  log  Qb  flow  estimated  by  regression  of  base  flow  measurements 
Sb^  &  Sr'    are  variances  in  square  log  units  of  estimates  from  true  value  of  low  flow 
characteristics 
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An  Infiltration  Index  Useful  in  Estimating  Low  Flow  Q-iaracteristics  of  Drainage  Basins 
(Armbruster,  J.T.,  1976a) 

This  investigation  concluded  that  the  basin  geology  was  an  overriding  factor  in  describing 
low  flow  characteristics.  Using  an  SCS  soil  type  based  infiltration  index,  the  author  of  this 
study  was  able  to  improve  standard  errors  in  existing  prediction  equations.  The  soU  index 
was  determined  by  percentage  of  basin  of  soil  type  A,  B,  C  or  D  and  applying  a  weighting 
of  10,  5,  1,  1  respectively  and  summing  for  the  entire  basin.  The  drainage  area  and  a 
precipitation  index  were  also  found  to  be  significant  factors  in  estimating  low  flows.  The 
final  equation  takes  the  following  form: 

log  7Qio  =  -4.658  +  1.068  log  DA  -h  2.183  log  MAP  +  .218  SI 

where  SI  is  Soil  Index 

Technical  Manual  for  Estimating  Low  Flow  Frequency  Characteristics  of  Streams  in  the 
Susquehanna  River  Basin  (Armbruster,  I.T.,  1976b) 

This  technical  report  adopted  the  ii\filti-ation  index  described  in  the  earlier  reference  by 
Armbruster  which  examined  step-by-step  procedures  for  low  flow  frequency  characteristics. 
In  addition,  distinctions  were  made  between  large  and  small  drainage  basins.  Large  basins 
have  drainage  areas  of  larger  than  5,120  km^  (2000  mi^).  The  parameters  used  in  the 
regressions  were;  draiiuge  area,  mean  annual  precipitation,  basin  surface  storage,  index  of 
relative  infiltration.  As  well,  an  index  coefficient  of  the  1  day  flow  to  the  3  day  flow  of  .93 
was  determined.  For  large  basins,  a  graph  of  correction  factors  for  channel  length  was 
produced  and  applied  to  the  regression  equation  as: 

log  (7Q20)  =  -2.029  +  1.052  log  DA  +  .753  log  (MAP)  +  CF 

where  CF  is  the  correction  factor  from  the  graph. 
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Estimates  of  Low  Flows  Using  Watershed  Climatic  Parameters  (Boyer,  P.G.  and  M.  Chang, 

1977) 

A  regional  analysis  was  vmdertaken  for  12  unregulated  mountainous  watersheds  with 
drainage  areas  ranging  from  164  km^  to  2345  km^  (64-916  mi^)  located  in  West  Virginia. 
There  was  considerable  effort  put  forward  in  developing  physical  watershed  parameters  for 
which  the  evap>orative  power  and  drought  pxjtential  were  related  logically.  Eighteen 
watershed  parameters  were  tested;  drainage  area,  length  of  perimeter,  mean  elevation,  mean 
latitude,  maximum  relief,  relief  ratio,  basin  slope,  form  factor,  main  channel  slope,  main 
channel  length,  drainage  density,  stream  frequency,  stream  order,  average  angle  of  tributaries, 
aspect,  length-width  ratio,  coefficient  of  compactness  and  percent  of  forest  cover.  Climatic 
parameters  used  for  analysis  were  7  day  mean  maximum  temperature  and  maximum 
consecutive  days  without  precipitation  for  seasonal  periods  of  August  to  October. 

The  final  model  used  3  watershed  and  2  climatic  parameters,  namely;  watershed  perimeter 
(CO),  the  main  channel  length  in  miles  (LEN),  the  Watershed  Form  Factor  (WF) 
(area/maximum  watershed  length^).  The  7  day  10  year  maximum  temperature  (T]o),  a 
fimction  of  mean  elevation  and  latitude,  and  the  September  10  year  maximum  consecutive 
rainless  days  (Rjo),  and  finally,  a  function  of  watershed  mean  elevation.  A  prediction 
equation  was  develop)ed  for  the  1:10  year  7  day  low  flow  with  a  standard  error  of  estimate 
of  about  30%  of  the  observed  mean: 

In  (7Q,o)  =      37.0469  +  0.03886  (CO)  -  0.0401  (LNTH)  -  3.6743  (WF)  -  .5771  Tjo  +  1.1665  Rjo 

It  was  concluded  that  use  of  climatic  parameters  which  index  drought  conditions  can  improve 
low  flow  estimates,  and  that  meaningful  low  flow  estimates  can  be  obtained  from  climatic 
and  watershed  parameters  or  watershed  parameters  alone  in  a  mountainous  humid  region. 

Analysis  of  Low  Stream  Flows  on  Cape  Breton  Island,  Nova  Scotia  (Environment  Canada, 
1978) 

This  investigation  included  frequency  analysis  using  Gumbel  IE,  for  1,  3,  7, 10,  30  day  low 
flow  durations  based  on  11  hydrometric  stations  with  drainage  areas  ranging  from  216  km^ 
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to  3.6  km^  (142  -  2.4  mi^).  Stations  with  periods  of  record  from  10  to  59  years  were  used  aitd 
those  with  man-made  regulation  were  deleted. 

It  was  determined  that  extreme  low  flow  estimates  for  the  short-term  hydrometric  statioi^s 
may  be  high.  Seasonal  analysis  found  higher  low  flows  for  the  summer  period.  It  was 
suggested  that  this  was  due  to  the  dimatic  and  physical  location  of  the  stations,  and  it  was 
noted  that  evaporative  power  may  also  be  important.  As  well,  a  physiological  data  base  was 
prepared  for  a  future  regression  study  of  13  parameters  relating  to  the  stations  examined  (see 
Table  2.1). 


Journal  of  the  Hvdraulic  Division  -  Characteristics  of  Low  Flows  (Task  Committee,  Surface 
Water  Hydrology,  1980) 

This  papjer  summarizes  the  findings  of  a  U.S.  committee  regarding  the  types  of  low  flow 
information  needed,  various  low  flow  characteristics  and  their  accuracy,  and  suggests  further 
analysis  and  data  collection  needs.  The  main  topics  were  discussed  according  to  gauged  sites 
and  vingauged  sites. 

a)  Gauged  Site 

The  most  widely-used  index  of  low  flow  in  the  U.S.  is  the  7  day  10  year  low  flow  (7Q,o). 
However,  it  was  also  found  tiiat  other  durations  of  from  1-273  days  have  been  based  on  data 
records  of  15-20  years  or  more. 

b)  Ungauged  Sites 

It  was  found  that  low  flow  characteristics  depend  largely  on  the  geology  of  the  basin  and  on 
losses  from  evapotranspiration.  For  example,  the  incorporation  of  basin  rock  type  into 
regional  analysis  of  low  flows  increased  prediction  accuracy  for  equations  developed  in 
Virginia.  A  regional  low  flow  forecasting  model  developed  (by  Wright,  1970)  for  southeast 
England  relates  lowest  mean  flow  for  the  forecast  year  to;  drainage  area,  mean  annual  runoff, 
geology  index,  mean  summer  catchment  rainfall  in  the  same  year  as  forecast  and  mean  winter 
catchment  rainfall  (October-March)  preceding  forecast  low  flow. 


2-7 
T  nw  Flow  Studies  (Institute  of  Hydrology,  1980) 

The  Institute  of  Hydrology  (U.K.)  analysed  low  flows  using  regression  techniques.  An 
innovative  feature  was  the  development  of  a  geologic  index  called  the  baseflow  index.  This 
index  can  be  determined  from  streamflow  measurements  (for  equation  see  Condie,  1986).  In 
ungauged  areas,  comparative  areal  methods  can  be  used  to  estimate  the  BH  or  it  can  be 
estimated  from  underlying  geology. 

Low  flow  duration  curves  were  regressed  ft-om  basin  parameters.  Then  a  comparative 
analysis  of  the  change  of  slope  for  different  n-day  durations  allowed  for  a  prediction  of  any 
n-day  low  flow.  In  the  study,  much  work  was  done  to  provide  non-dimensional  graphs  and 
to  linearize  them  for  simplifying  the  estimation  procedures. 

The  regional  equatioi\s  tcike  the  form: 

Q,5(10)    =QBn  +  QMAP-C3 
andQ«(10)    =  C4  BFI  +  C5  LNTH  -  Q 

where    MAP  =  mean  annual  precipitation 

LNTH  =  stream  length  for  region  five. 

As  well,  the  slope  is  regressed  from  parameters  as: 

logio  GRADMEN  0,5  =  07  MAP  -  Cg  Q95(10)  -  Q 

These  two  regression  equations  were  then  used  to  estimate  Q„  (n)  for  any  n<iay  duration  by 
the  formula: 

Q,5  (n)  =  Q,5  (10)  [1  +  (D-10)  GRADMEN  Q^] 

Similar  relationships  were  developed  for  mean  annual  minimum  10-day  low  flow  MAM  (10). 
Included  in  the  report  describing  this  five  year  study  is  a  fairly  extensive  Uterature  review 
which  is  paraphrased  below. 
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Review  of  the  Institute  of  Hydrology  Background  Review  flnstitute  of  Hydrology,  1980) 

"Many  examples  worldwide  employ  catchment  characteristics  indexing  catchment  size,  shape 
and  climate  and  most  refer  to  the  additional  need  to  consider  the  geology  of  the  catchment 
in  order  to  explain  more  of  the  Vcuiability  of  low  flows.  Wright  (1970)  presents  a  prediction 
equation  for  the  lowest  flow  in  a  year  using  catchment  slope  and  area  as  independent 
variables.  Riggs  (1973),  Wesche  et  al  (1973),  and  Skelton  (1974),  make  use  of  the  cross- 
sectional  properties  of  streams  to  correlate  low  flows  with  neighbouring  gauged  catchment 
data.  For  the  United  Kingdom,  general  information  on  behaviour  of  British  rivers  is  given 
by  Ward  (1968)  and  also  Rodda  (1976).  The  interaction  between  surface  and  groimdwater 
for  major  aquifers  is  discussed  in  Ineson  et  al  (1965).  Much  design  work  has  been  based  on 
the  lowest  recorded  flow.  Flow  duration  curves  and  their  application  in  the  U.K.  are  given 
in  Hoyle  (1963).  They  are  discussed  by  Boulton  (1965)  where  the  concept  of  a  Minimum 
Acceptable  Row  is  described  and  the  factors  likely  to  enter  into  its  calculations  are 
considered." 

Flood  Prone  Area  and  Low  Flow  Analysis  Study  in  the  Kawartha  Region  (MacLaren 
Plansearch,  1981) 

This  study  developed  a  regression  equation  for  predicting  low  flow  estimates  in  the 
Kawarthas  region.  The  regression  equations  took  the  form: 

-^  =  1.64  X  I0-'  SP"'  ALS^"« 
DA 


-^  =  6.90  X  10^  S/'^  ALS-^-^'^ 
DA 


where         qi  =  average  minimum  daily  flow 

Qs  =  average  minimum  summer  flow 

DA  =  drainage  area 

SI  =  index  of  average  soil  moisture  holding  capacity 

ALS  =  %  area  covered  by  swamp  or  wetland 
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Application  of  Statistical  Low  Flow  Analysis  as  a  Basis  for  Water  Quality  Planning  (Rubach, 
H.,  1982) 

Two  practical  examples  for  applications  of  statistical  low  flow  analysis  as  related  to  water 
quality  management  were  discussed.  In  the  first  example,  water  quality  parameters  were 
described  as  a  function  of  runoff  of  defined  frequency.  This  was  used  as  information  to 
determine  the  necessary  amount  of  low  flow  augmentation.  In  the  second  example,  the 
necessary  size  of  storage  volumes  desired  for  low  flow  augmentation  was  estimated.  The 
storage  volumes  depended  on  the  desired  minimal  runoff  vcilue,  the  frequency  and  the 
duration  of  the  augmentation  process.  To  generate  information  about  mean  low  flows  for 
defined  durations  and  frequencies,  a  duration-runoff-frequency  fimction  was  determined.  On 
the  basis  of  daily  runoff  values,  the  mean  low  flow  was  calculated  for  a  set  of  different 
durations  of  low  flow  pjeriods.  For  each  duration  step,  the  yearly  lowest  value  of  mean  low 
flow  in  the  period  of  this  duration  was  used  as  an  element  of  a  series  of  extreme  values.  A 
frequency  distribution  was  fitted  to  this  empirical  series  of  extreme  values  to  extrapolate  low 
flows  of  defined  frequencies. 

Low  Flow  Frequency  Analysis  (Program  LOFLOW)  (Condie,  R.  and  L.  Cheng,  1983) 

The  Gumbel  Type  III  distribution  is  discussed  as  well  as  conceptual  problems  found  from  the 
Condie,  Nix  1975  study.  It  was  found  that  the  Gumbel  in  distribution  becomes  unstable  if 
the  sample  has  a  skewness  of  less  than  -1.08.  This  study  concluded  that  the  Gumbel  III 
distribution  should  be  used  as  the  basic  method  for  low  flow  frequency  analysis.  As  an 
alten\ative  for  samples  with  large  negative  skewness,  it  was  suggested  that  a  3-parameter  log 
normal  distribution  should  be  used. 

The  only  criteria  of  successful  fit  to  the  Gumbel  III  is  that  the  lower  boundary  parameter 
should  not  be  larger  than  the  smallest  member  of  the  sample.  Parameter  estimation  should 
proceed  in  the  order  of;  maximum  likelihood,  smallest  observed  drought  and  by  moments. 
The  program  LOFLOW  was  written  to  analyse  samples  of  low  flows  taking  these  concerns 
into  account. 
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Computation  of  the  Base  Flow  Index  program  (BFINDEX)  (Swan,  P.  and  R.  Condie,  1983) 

It  was  determined  that  low  flows  are  very  dependent  on  catchment  geology,  and  the  Base 
Flow  Index  was  suggested  as  a  practical  technique  to  provide  an  index  of  this  parameter.  A 
program  was  written  to  Ceilculate  a  base  flow  index  (BFI)  as: 

BFI  =   total  volume  of  base  flow/yr 
total  volume  of  nmoff /yr 

The  required  data  was  obtained  from  records  published  by  the  Water  Survey  of  Canada,  and 
averaged  BFI  values  were  computed  over  a  period  of  several  years. 

The  main  problem  with  this  technique  was  determined  to  be  the  requirement  for  discharge 
measurements.  However,  it  was  pointed  out  that  the  BFI  parameter  was  found  to  be  fairly 
stable  over  both  wet  aiui  dry  years  and,  it  was,  therefore,  concluded  that  it  could  be 
estimated  from  as  little  as  one  to  two  years  of  discharge  measurements. 

Effects  of  Man  on  Low  Flows  (Riggs,  H.C.,  1984) 

This  investigation  examined  the  climatic,  geologic  and  topxjgraphic  influences  on  low  flow 
characteristics  in  view  of  man's  influence.  Man-made  changes  that  affect  low  flows  were 
identified  as  changes  in  vegetation,  drainage  practices,  diversions,  return  flows,  groundwater 
pumping,  sewage  treatment  plants  and  other  effluent  discharge. 

It  was  pointed  out  that  the  aquifer  characteristics  are  also  important  to  low  flows.  It  was 
postulated  that  evapotranspiration  and  changes  to  ground  cover  can  likely  be  related  to  low 
flow  increases.  Urbanization  is  not  a  significant  factor  since  pluses  (i.e.  sewer  leakage)  and 
minuses  (i.e.  reduced  infiltration)  tended  to  even  out.  Reservoir  regulation  usually  increases 
low  flows  due  to  minimum  outflow  requirements.  Diversions  are  highly  basin  dependent 
and  are  diiflcult  to  assess. 

It  was  noted  that  swamps  apparently  had  very  smaU  low  flows  in  extreme  drought  cases. 
Therefore,  it  was  postulated  that  changes  to  swamps  (e.g.  by  infilling,  etc.)  might  not  have 
a  significant  impact  on  downstream  low  flow  characteristics.  Pumping  can  have  a  large  effect 
depending  on  duration  location  and  basin  geology.  The  three  main  parameters  for  low  flow 
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are  related  by;  climate  providing  input  and  output,  geology  and  topography  influendng 
infiltration  to  aquifers,  and  then  geologic  parameters  affecting  transport  to  streams. 

Man's  influence  on  precipitation,  temp)erature  and  irrigation  are  not  likely  to  affect  low  flows 
while  changing  vegetation  cover  (evapotranspiration),  regulation  and  diversion  are  likely  to 
have  the  most  evident  effects  on  low  flow. 

Hydrologic  Oesign  Methodologies  for  SmaU  Scale  Hydro  at  Ungauged  Sites  (Acres,  1985) 

This  study  develojjed  a  regression  relationship  between  basin  physiographic  parameters  and 
"turbinable  flow"  which  was  related  to  exceedance  curves.  Exceedance  curves  were  non- 
dimensionalized  by  dividing  flow  values  by  the  mean  annual  flow.  For  ungauged  sites,  a 
regression  equation  provided  estimates  of  mean  annual  flow  from  values  of  MAR  and  area. 

MAF  =  0.0000369  MAR-DA 

This  was  then  related  to  a  regional  indexed  curve  to  determine  values  of  exceedance  by 
examining  the  regional  mean  annual  flow  to  the  sight  mean  annual  flow. 

As  weU,  two  other  methods  were  compared,  proration  and  ARMA  (Autoregressive  Moving 
Average  Methodology).  It  was  determined  that  proration  provided  the  best  practical  method 
for  determining  discharge.  The  pn-oration  formula  took  the  form: 

Qjk  =  Qlijk  (DA  MAR)  /  (DAI  MARI) 
where  i  =  year,  j  =  month,  k  =  day 

The  ARMA  method  was  used  to  synthesize  daily  and  monthly  flows  from  required 
parameters  within  the  basin  over  a  time  series. 

Median  Drought  Flows  at  Ungauged  Sites  in  Southern  Ontario  (Pilon,  P.  and  R.  Condie,  1986) 

A  methodology  was  suggested  and  explored  for  predicting  median  drought  flows  at 
ungauged  sites  in  Southern  Ontario.  Minimum  average  flows  for  durations  of  3, 10, 14,  21 
and  28  consecutive  days  were  determined  each  year  in  the  time  period  June  1  to  May  31  for 


2-12 

74  stations  in  Southern  Ontario.  A  regression  equation  was  developed  for  predicting  the 
median  10  day  low  flow  using  drainage  area  (DA),  mean  annual  precipitation  (MAP),  and 
base  flow  index  (BFI).  The  precipitation  parameter  was  not  significant  and,  therefore,  not 
included  in  the  equation.  The  standard  error  of  the  low  flows  was  found  to  be  37%  and  the 
equation  was  in  the  form: 

MED(IO)  =  0.000075  DA  +  1.515  BFI  -  0.356 
The  median  10  day  low  flow  was  then  used  to  produce  a  slop>e  relationship  to  relate  other 
low  flow  durations  to  the  MED(IO)  value.  This  equation  is  expressed  as: 

GRADMEN  =  2.632  MED(IO)  +  0.46399 
Where  GRADMEN  is  the  slope  of  the  median  flow  duration  relationship. 

A  problem  noted  by  the  authors  is  the  detennination  of  BFI  for  ungauged  sites.  The 
parameter  is  used  as  an  index  of  catchment  geology. 

Drought  Flows  and  Receiving  Water  Assessment  (Logan,  L.,  1986) 

This  study  examined  low  flows  in  an  environmental  context  v^th  respect  to  waste  water  and 
failures  to  satisfy  water  quality  requirements.  Low  flow  design  estimates  for  Central  - 
Southeastern  Ontario  were  reported  to  be: 


Central 

Southeastern 

(mVs/km^ 

mVs/km^ 

7Qa, 

1.0 

0.4           ) 

7Q,o 

1.1 

0.6          )     ±  .3  mVs/km" 

7Q, 

1.8 

1.7          ) 

The  MOE  guideline  on  drought  flow  design  is  7Q20.  This  implies  that  lower  low  flows  will 
cause  water  quality  violation  in  BOD  and  other  wastes. 

Approach  for  Frequencv  Analysis  of  Multiyear  Drought  Durations  (Lee,  K.S.,  J.  Sadeghipour, 
and  J. A.  Dracup,  1986) 

An  approach  for  frequency  analysis  of  multiyear  drought  durations  was  presented.    The 
problem  of  sample  size  is  discussed  with  a  technique  for  smoothing  frequency-curve 


2-13 

irregularities  of  drought  durations.  Piecewise  linear  regressions  for  the  linearized  termination 
probability  curve  define  the  parameters  of  the  logistic  model.  The  model  then  can  be  used 
to  determine  the  return  periods  of  multilinear  droughts.  This  report,  although  interesting, 
is  not  directly  applicable  to  the  present  investigations. 

Regional  Flood  Frequency  Analysis  for  Ontario  Streams,  Volume  2,  Multiple  Regression 
Method  (Moin,  S.M.A.,  and  M.A.  Shaw,  1986 

This  report  primarily  concerns  prediction  of  flood  peaks  although  it  includes  a  discussion  of 
the  Base  Flow  Index  which  may  be  relevant  for  the  present  investigations.  The  BFI  was 
calculated  and  mapped  for  2  flood  regions  identified  in  the  Province. 

The  regression  equations  wdth  BFI  for  predicting  flood  flows  provided  predictions  superior 
to  those  derived  using  the  secondary  regression  equatioris  v^thout  the  BFI  parameter.  The 
latter  was  found  to  be  second  only  to  drainage  area  as  a  significant  parameter  in  describing 
vciriance  of  flood  flows. 

The  BFI  was  foimd  to  be  extremely  sensitive  to  basin  storage  conditions  and,  therefore,  could 
not  be  used  in  the  primary  regression  equations  for  the  northern  region  due  to  the  large 
nimiber  of  lakes  and  swamp)s. 

Correlations  with  SCS  soil  groups  indicated  that  BFI  could  be  approximated  as  follows: 

Soil  Group  A  >  0.60± 

Soil  Group  B  0.45± 

Soil  Group  C  0.30± 

Soil  Group  D  <  0.25± 

Low  Flow  Studies  in  Scotland  (Institute  of  Hydrology  and  Department  of  the  Environment, 
1986) 

The  report  described  a  method  of  estimating  monthly  and  seasonal  flow  diiration  curves  at 
sites  where  litfle  or  no  flow  data  are  available.  The  method  was  developed  to  estimate  the 
95  percentile,  10  day  flow  from  the  aimual  flow  duration  curve,  Q^  (10),  10  day  average  flow 
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Low  Flow  Studies  in  Scotland  flnstitute  of  Hydrology  and  Department  of  the  Environment, 
1986) 

The  report  described  a  method  of  estimating  monthly  and  seasonal  flow  duration  curves  at 
sites  where  Uttle  or  no  flow  data  are  available.  The  method  was  developed  to  estimate  the 
95  percentile,  10  day  flow  from  the  annual  flow  duration  curve,  Q,^  (10),  10  day  average  flow 
exceeded  by  95%  of  10  day  average  discharge,  Qjj  (10)  was  estimated  in  units  of  percentage 
of  average  flow  in  m^/s  by  using  regression  equations  developed  for  different  regions  using 
the  Base  Flow  Index  (BFI)  and  standard  period  (1941-70)  armual  average  rainfall  (SAAR)  or 
defined  as  MAP  (Mean  Annual  Precipitation)  in  the  present  study.  The  generalized  form  of 
the  regression  equation  was: 


For  Flow  estimation  at  the  ungauged  sites,  regression  equatior\s  were  revised  for  each  of  these 
regions  for  estimating  Q95  (10)  and  MAM  (10)  (mean  annual  10  day  minimum)  from  the  Base 
Flow  Index  (BFI),  Standard  Annual  Average  Rainfall  (MAP),  and  Proportion  of  Catchment 
Covered  by  a  Lake  or  Reservoir  (FALAKE)  or  ACLS  (Area  controlled  by  lakes  and  swamps 
as  defined  in  the  present  study).  The  generalized  low  flow  regression  equations  are  given 
below: 

^/Q95aO)  =  flo  +  Oi  y/BFI  +  Oj  y/MAP  -  a,  y/XCU 
v/A£AM(10)  =  a^  *  a^  y/SFl  -  a^  y/ACLS 

The  negative  regression  coefficient  of  ACLS  does  not  imply  that  catchments  with  lakes  have 
lower  low  flows  than  those  without.  This  apparent  paradox  is  resolved  by  recalling  that  the 
attenuating  effect  of  a  lake  on  the  downstream  hydrograph  will  greaUy  increase  the  BFI.  This 
results  in  higher  BFI  in  catchments  with  lakes  for  a  given  (295(10)  than  in  lake-free 
catchments  and  so  this  increased  BFI  is  compensated  in  the  equation  by  a  negative  coefficient 
of  ACLS.  The  same  phenomenon  was  noted  in  the  equation  developed  for  predicting 
MAM(IO).  The  other  notable  feature  of  the  equatior^  is  that  SAAR  is  a  useful  explanatory 
variable  for  (295(10)  but  not  for  MAM(IO).  This  was  attributed  to  the  finding  that  (295(10) 
was  higher  in  wet  than  in  dry  areas  having  the  same  BFI  but  that  MAM(IO)  was  independent 
of  rainfall  in  most  regions. 
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A  Comparison  of  Methods  for  Estimating  Low  Flow  Characteristics  of  Streams  (Tasker,  CD., 
1987) 

Four  methods  for  estimating  the  7  day,  10  year  and  7  day,  20  year  low  flows  for  streams  were 
compared  by  the  bootstrap  method  (Log-Pearson  III,  Weibull,  Box-Cox  and  Log-Boughton). 
The  bootstrap  method  is  a  Monte  Carlo  techruque  in  which  random  samples  are  drawn  from 
an  imspediied  sampling  distribution  defined  from  observed  data.  The  non-parametric  nature 
of  the  bootstrap  makes  it  suitable  for  comparing  methods  based  on  a  flow  series  for  which 
the  true  distribution  is  uiUcnown.  It  was  found  that  when  all  the  data  were  used,  all  four 
methods  gave  reasonably  dose  estimates.  However,  in  general,  the  methods  based  on 
hypothetical  distributions  (the  LP  111  and  Weibull  methods)  performed  better  in  terms  of 
mean  square  error  than  did  the  Box-Cox  transform  method  or  the  Log-Boughton  method 
which  is  based  upon  a  fit  of  plotting  positions.  The  LP  in  method  performed  marginally 
better  than  the  Weibull  except  for  very  short  record  lengths  where  both  methods  were  found 
to  be  about  equal.  The  authors  concluded  that  the  results  reported  did  not  extend  over  the 
entire  range  of  low  flows  and  did  not  imply  that  one  particular  distribution  was  better  than 
another  in  describing  the  distribution  of  low  flows. 

Low  Flow  Characteristics  in  Ontario  (Cumming  Cockbum  Limited,  1989) 

The  study  was  undertaken  to  update  the  statistical  low  flow  analyses  for  all  the  five 
administrative  regions  and  to  produce  extreme  value  analyses  for  suitable  stream  gauging 
locations  (1,  3,  7,  15,  30  day  durations)  with  record  lengths  greater  than  10  years.  The 
minimum  aimual  consecutive  7  day  average  low  flow  was  determined  for  each  year  and  the 
corresponding  set  of  consecutive  7  day  average  low  flows  represented  an  extreme  value  series 
to  which  a  theoretical  extreme  value  distribution  was  fit.  It  was  noted  that  previous 
investigations  generally  utilized  the  Gumbel  in  distribution.  In  addition,  statistical 
characteristics  of  average  low  flow  values  for  all  of  Ontario  were  compared  to  those  for  each 
region.  It  was  noted  that  a  higher  than  expected  number  of  station's  data  illustrate  some 
trend  and  persistence. 
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Environmental  Ontario  Research  Project  :  Assessment  of  the  Biologically  Based  Low  Flow 
Analyses  Technique  (Gumming  Cockbum  Limited,  1989) 

This  study  analysed  and  compared  the  existing  hydrologic  methodology  for  defining  low 
flows  (extreme  value  analysis)  to  a  biologically  based  methodology  developed  by  the  U.S. 
EPA  (Environmental  Protection  Agency).  The  report  concluded  that  the  biological  low  flow 
model  could  successfully  be  applied  to  river  systems  in  Ontario.  The  ZQzo  hydrologically 
defined  low  flow  is  sunilar  to  the  chronic  design  flow  in  terms  of  flow  magnitude  and  the 
exp)ected  number  of  excursions  from  water  quality  objectives.  The  report  recommended  that 
the  MOE  should  continue  use  of  the  7Q20  as  the  design  low  flows  for  the  Province  of  Ontario. 

Generalized  Low  Flow  Frequency  Relationships  for  Ungauged  Sites  in  Massachusetts  (Vogel, 
R.M.  and  G.N.  KroU,  1990) 

Generalized  regional  regression  equations  for  estimating  the  n-day,  y-year  low  flow 
discharge,  nQy,  at  ungauged  sites  in  Massachusetts  were  developed  for  n  =  3,  7, 14  and  30 
days.  A  2-parameter  log  normal  distribution  was  fit  to  sequences  of  annual  minimum  n-day 
low  flows  and  the  estimated  parameters  of  the  log  normal  distribution  were  then  related  to 
two  drainage  basin  characteristics;  drainage  area  and  relief  (for  23  sites  with  natural  flow). 
It  was  pointed  out  that  many  previous  investigations  had  developed  regional  models  for 
estimating  low  flow  statistics  at  imgauged  sites  from  readUy  available  geomorphic,  geologic, 
climatic  and  topographic  characteristics  (e.g.  Thomas  and  Benson,  1990).  Usually  these 
models  take  the  form: 

nQy  =  b„X^'^X^^^)^^3 

where  nQy  is  obtained  from  gauged  streamflow  records,  the  X;  are  easily  measiored  drainage 
basin  and  climatic  characteristics,  and  the  b;  are  parameter  estimates  obtained  from  the  use 
of  multivariate  regression  procedures.  Vogel,  et  al  (1989)  presented  a  regional  low  flow 
model  that  describes  the  time  response  of  runoff  from  a  watershed  during  dry  periods. 
Describing  the  watershed  as  a  system  of  linearized  Dupuit-Boussinesg  aquifers,  they  found 
that  watershed  runoff  (Q)  at  time  t^  after  the  hydrograph  peak  may  be  described  by  the 
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Q,  =  4k  (DA)  (Hd)^  Kb* 

where  k  is  the  hydraulic  conductivity  of  the  aquifer,  DA  is  drainage  area,  H  is  watershed 
relief,  d  is  drainage  density,  K^  is  the  dimensionless  baseflow  recession  constant  and  t  is  time. 
The  product  SLP  is  an  approximation  to  the  average  basin  slope.  The  basin  relief,  H,  is 
defined  as  the  difference  in  elevation  between  the  basin  summit  and  the  basin  outlet  where 
the  basin  sununit  elevation  is  the  average  of  the  highest  peak  and  the  two  adjacent  peaks  on 
either  side  of  it.  Unfortunately,  accurate  estimates  of  k  and  kj,  were  fovmd  to  be  difficult  to 
obtain  at  ungauged  sites;  hence  those  characteristics  were  ignored.  The  inclusion  of  drainage 
density,  d,  did  not  result  in  improvements  of  estimates  of  low  flow  characteristics  at 
imgauged  sites,  hence  that  basin  characteristic  was  ignored,  hence  this  study  (Vogel  and 
Knoll,  1990),  assumed  that  low  flows  could  be  characterized  by  watershed  area  and  basin 
relief  alone.  The  low  flow  predictions  from  this  model  were  comparable  with  previous 
regional  low  flow  investigations  in  terms  of  average  prediction  errors.  It  was  emphasized 
that  further  reductions  in  the  prediction  errors  associated  with  such  regional  low  flow  models 
will  likely  result  from  the  inclusion  of  additional  watershed  characteristics  such  as  basin 
hydraulic  conductivity,  average  basin  slopje  and  the  baseflow  recession  constant.  The  study 
also  compared  the  use  of  generalized  least  squares  and  (OLS)  regression  procedures.  The 
study  concluded  that  the  GLS  and  OLS  regression  procedures  led  to  almost  identical  regional 
regression  model  parameter  estimates  although  GLS  procedures  led  to  only  marginal  gains 
in  the  prediction  errors  associated  with  low  flow  regional  regression  equations. 

Regional  Analysis  of  Low  Flow  Characteristics  (Southwestern  and  West  Central  Regions) 
(Cumming  Cockbum  Limited,  1990) 

The  study  develop>ed  regional  regression  equations  based  on  Drainage  Area  (DA),  Base  Flow 
Index  (BFI),  stream  length  (LNTH),  other  parameters,  for  example.  Mean  Annual  Runoff 
(MAR)  and  Mean  Annual  Snowfall  (MAS)  were  also  examined.  The  generalized  regression 
equation  developed  was  of  the  following  form: 

T^Qy  =  ao  +  ai  DA^  +  ai  BFl^  +  a^  LNTH^ 

where  nQy  is  the  n-day  average  low  flow  occurring  once  in  y  years. 
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The  study  concluded  that  the  regression  utilizing  E>rainage  Area,  Base  Flow  Index  and  Stream 
Length  as  independent  parameters  provide  reasonable  estimates  of  low  flow  values.  It  was 
also  found  that  the  equations  for  estimating  low  flows  were  sensitive  to  variations  in  values 
of  Base  Flow  Index  (BFI)  and  Length  (LNTH).  It  was,  therefore,  concluded  that  care  should 
be  taken  in  estimating  these  parameters  for  ungauged  watersheds.  For  estimation  of  the  BFI, 
it  was  concluded  that  available  maps  of  BFI  isolines  could  be  used  or  short  term  monitoring 
(2  years)  could  provide  good  estimates.  Finally,  it  was  also  concluded  that  the  Index  or 
graphical  and  regression  methods  provided  superior  estimates  of  low  flows  compared  to 
those  obtained  solely  from  station  proration  of  imit  area  low  flows  from  gauged  sites. 

Summary  of  Literature  Review 

In  summary,  it  was  foimd  that  many  investigators  have  analysed  low  flow  characteristics 
over  durations  from  1  to  273  days  in  length,  for  return  periods  of  1  to  100  years.  In  general, 
it  has  been  determined  that  the  Gumbel  IE  and  Log-Pearson  HI  distributions  best  fit  the 
samples  of  low  flow  data  for  gauged  streams  (Tasker,  1987,  Condie  and  Nix,  1975,  Matalas, 
1965  and  Table  2.1).  Several  investigations  considered  techniques  for  predicting  selected  low 
flow  characteristics  while  very  few  investigations  have  considered  techniques  for  developing 
regiorulized  flow  duration  curves. 

Parameters  foimd  to  be  significant  in  regression  equatioris  generally  tended  to  be  drainage 
area,  base  flow  index,  mean  annual  precipitation,  area  controlled  by  lakes  and  swamps, 
watershed  relief,  stream  length,  mean  annual  snowfall,  mean  annual  evaporation, 
groundwater  fluctuations,  mean  annual  runoff  and  soil  index.  However,  some  parameters 
that  were  not  found  to  be  highly  correlated  in  one  region  were  found  to  be  important  in  other 
regions  (Institute  of  Hydrology,  1980). 

With  reference  to  the  available  literature,  it  was  decided  to  include  the  follovsdng  parameters 
in  development  of  the  physiographic,  meteorologic  and  hydrometric  data  base;  drainage  area 
(DA),  mean  emnual  snow^cill  (MAS),  mean  annual  precipitation  (MAP),  mean  annual  runoff 
(MAR),  mean  annual  evaporation  (EVA),  area  controlled  by  lakes  and  swamps  (ACLS),  slope 
(SLP),  stream  length  (LNTH),  Base  Flow  Index  (BFI)  and  maximum  groundwater  fluctuations 
(GW). 


TABLE  2.1 
SUMMARY  OF  LITERATURE  REVIEW 


Study 

No.  of 
Stations 

Distri- 
bution 

Method 

Parameters 
Examined 

General  Form  of  Equations 

Wright,  1970 

Regression 

DA,BS,MAP 

MAM  =  f(DA,  BS) 
MAM  =  f{DA,  BS,  MAP) 

Osbom,  1974 

20 

Index 
Regression 

Q7L1P,P,DA,H,CL,DD 

log  7Qj  =  f(DA,  H,  Q7LIP,  P) 

Taster,  1975 

LNHI 

Regression 

7Q,o  =  KDA,  MAP,  T) 

Armbruster,  1976a 
Annbruster,  1976b 

104 

Regression 

DA,CS,BL,EL,FCMAP, 
ACLS,MAS,SI,MAP,CF,DA 

log  7Q,o  =  f(DA,  MAP,  SI) 
7Q,o  =  KDA,  MAP,  CF) 

Boyer,  1977 

12 

G 

Regression 

DA,CO^L,LAT,LNTH,WF, 
T,o,R,o 

In  (7Q,o)  =  f(CO,  LNTH,  WF,  T,o,  R,o) 

Env.  Can.,  1978 

11 

G 

DA,R,AS,AL,ALS,OB,LEN, 
DD,SLP,EL,MAP 

Correlation  Matrix  only  produced 

Inst,  of 
Hydrology,  1980 

G 

Regression 

Q„(10),Bn,DA,MAP,LNTH 

Q^dO)  =  F((Bn)^,  (DA)^,  (MAP)-^) 
Qsb(IO)  =  F((Bn)-',  (DA)-=,  (LNTH)^) 
(lOOj)  =  fCBH) 
(lOOj)  =  fCBH,  MAP) 

MacLaren 
Plansearch,  1981 

9 

Eye 
Fit 

Regression 

DA,SI,ALS 

Q  =  KDA,  SI,  ALS) 
Q.  =  KDA,  SI,  ALS) 

Condie,  1983 

74 

G 

Regression 

DA,MAP,Bn 

MED(IO)  =  f(DA,  BFD 
GRADMEN  =  f  ((MED(10))^) 

Inst,  of 
Hydrology,  1986 

232 

LN 

Regression 

Bn,MAP,ACLS 

(095(10)**  =  KBFl,  MAP,  ACLS) 
(Q95(10)^  =  KBFl,  ACLS) 

Vogel  &  KroU, 
1990 

23 

LN 

Regression 

DA,H 

7Qy  =  f(DA,  H) 

Cumming 
Cockbum,  1990 

65 

3PLN 

Regression 

DA,Bn,LNTH 

nQy  =  f  (DA,  BH,  LNTH) 

2-19 


Recent  investigations  (Gumming  Cockbum  Limited,  1990)  also  recommended  that  the 
Ministry  of  the  Environment  should  continue  use  of  the  7Qjo  as  a  prime  indicator  of  low 
flows  for  the  Province  of  Ontario. 
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3.0  DATA  BASE 

3.1  General 

The  Central  and  Southeastern  regions  identified  by  the  Ministry  of  Environment  define  the 
study  area.  One  hundred  and  twenty-five  hydrometric  stations  were  identified  in  the  region 
for  this  analysis  based  upon  the  period  of  record  and  other  station  characteristics.  Low  flow 
characteristics  for  73  of  these  stations  were  analysed  to  develop  estimation  techniques. 
Twenty-one  of  those  remaining  were  set  aside  for  testing  of  results.  The  remaining  stations 
were  omitted  from  the  study  based  on  selective  screening  criteria.  (Station  screening  and 
selection  is  discussed  in  Section  3.2.2.) 

Low  flow  characteristics  at  each  station  were  retrieved  from  our  inhouse  computerized  data 
base.  The  methodology  for  analysing  extreme  low  flows  and  determining  relevant  recurrent 
intervals  is  discussed  in  the  Low  Flow  Characteristics  in  Ontario  report,  Cumming  Cockbum 
Limited,  Ministry  of  the  Environment,  1988  (Appendices  A,  B  and  C)  and  the  results 
summarized  on  Figures  3.1  and  3.2  of  this  report.  Relevant  physiographic  and  meteorological 
characteristics  were  also  determined  as  discussed  in  Section  33. 

3.2  Extreme  Value  Analysis 

Moving  average  low  flows  (n-day)  were  determined  and  extracted  for  each  year  of  the 
available  data  base.  Average  extreme  low  flows  were  determined  and  extracted  for  the  1,  3, 
7, 15  and  30  day  durations  and  are  available  as  part  of  the  background  files. 

An  extreme  value  aiialysis  was  then  tmdertaken  for  each  of  the  1,  3,  7,  15  and  30  day 
durations  for  each  of  the  stations.  The  results  of  the  analysis  are  presented  elsewhere 
(Cumming  Cockbum  Limited,  1988).  The  SPSS  program  was  then  used  to  produce  general 
statistics  of  the  data  base  including  the  mean,  standard  deviation  and  coefficient  of  skew  of 
the  available  low  flow  samples  for  different  durations.  These  general  statistics  are 
summarized  in  Tables  3.1  and  3.2. 
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3^.1      Summary  Maps  of  Gauged  Infonnation 

Selected  low  flow  c±iaracteristics  were  extracted  and  sunimarized  on  a  map  for  the  study  area 
(see  Figures  3.1  and  3.2).  The  station  locations  are  identified  by  a  A  located  on  the  rivers 
at  the  point  of  discharge  measurement  and  the  selected  data  is  summarized  in  an  information 
box. 

The  stations  (on  Figures  3.1  and  3.2)  are  identified  by  the  7  digit  Water  Survey  Nimiber  and 
this  is  followed  by  a  regulation  code.  The  "R"  indicates  that  data  collected  at  the  station  is 
affected  by  regulation,  "N"  means  the  station  data  are  natural  or  non-regulated.  The  boxes 
on  the  left  from  the  top  refer  to  the  average  7  day  flow  (m^/s)  with  a  recurrence  interval  of 
2,  5, 10  and  20  years,  followed  by  the  minimum  average  7  day  flow  and  the  period  of  record 
for  the  station.  The  values  shown  on  the  right  are  the  flows  (m^/s)  equalled  or  exceeded  for 
the  available  period  of  record  5, 50,  75,  95  and  99  percent  of  the  time  followed  by  the  station 
drainage  area  in  km^ 

Station  names  are  also  listed  along  with  the  station  numbers  for  identification  purposes  on 
Figures  3.1  and  3.2. 

322      Screening  Results 

The  stations  which  were  not  considered  to  possess  a  suitable  low  flow  data  bcise  for  the 
purposes  of  this  investigation  were  not  included  in  the  regionalization  analysis.  These 
stations  were  excluded  for  the  following  reasons: 

i)  Multiple  Stations 

It  was  determined  that  large  river  systems  may  have  several  streamflow  gauges  at  several 
locations  along  the  channel.  The  use  of  all  such  highly  correlated  data  could  adversely  bias 
the  development  of  regionalization  methods.  Therefore,  oiUy  representative  gauges  for  such 
streams  were  retained  for  further  analysis.  A  few  of  the  multiple  stations  were  used  for 
testing  the  developed  techniques. 
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ii)  Heavily  Regulated 

Low  flows  on  heavily  regiilated  streams  (i.e.  St.  Lawrence,  Ottawa  and  Trent  River  systems) 
are  affected  through  the  use  of  reservoirs  for  low  flow  augmentation.  Since  the  main  focus 
of  this  study  is  to  produce  a  method  for  investigating  low  flows  for  ungauged  watersheds, 
stations  indicating  a  high  degree  of  regulation  were  excluded  from  the  data  base.  The  Water 
Survey  of  Canada  identifies  gauges  as  regulated  (R)  or  natural  streams  (N).  The  degree  of 
regulation  is  not  quantified  for  these  regions.  Of  the  125  stations,  75  were  identified  as 
"regulated"  to  some  degree.  Therefore,'  objective  screening  of  the  "regvilated"  stations  was 
undertaken  (i.e.  stations  on  the  Trent,  Ottawa  and  St.  Lav^ence  system  were  removed).  The 
remaining  "regulated"  stations  were  retained  and  assigned  a  regulatory  code  representing 
regulation  or  non-regulation.  Subsequent  simple  correlation  analysis  (between  the  regulatory 
codes  and  the  low  flow  statistics)  confirmed  that  these  stations  could  be  retained  in  the  data 
set  due  to  insignificant  levels  of  correlation.  Future  studies  should  investigate  and  quantify 
degrees  of  regulation  and  p>ossible  techniques  to  deregulate  flow  series. 

iii)  Statistical  Tests 

Statistical  data  analysis  tests  were  undertaken  as  described  in  more  detail  in  Regional 
Analysis  of  Low  How  Characteristics  for  Southwestern  and  West  Central  Regions  by 
Cumming  Cockbum  Limited  (1990).  The  available  test  statistics  were  recently  made  available 
as  part  of  the  LFA  (PUon  and  Jackson,  1987)  low  flow  package. 

In  general,  it  was  found  that  a  significant  number  of  stations  "feiiled"  the  non-parametric  tests. 
Therefore,  taken  over  the  entire  data  base,  application  of  these  tests  has  indicated  that  the 
available  data  base  of  extreme  low  flows  may  exhibit  some  trend  and  dependence  with  some 
possibility  of  non-random  characteristics. 

The  data  were  further  analysed  by  subdivision  of  the  available  data  set  according  to  length 
of  record  (i.e.  >20  years  and  <20  years)  and  according  to  regulation  code.  However,  it  was 
foimd  that  neither  the  length  of  record,  nor  the  pxassible  effects  of  regulation  could  accoimt 
for  the  conclusions  of  the  test  results.  One  explanation  could  be  that  the  available  record 
lengths  are  too  short  to  permit  reasonable  application  and  interpretation  of  these  non- 
parametric  test  results.  Another  explanation  could  relate  to  seasonal  effects  on  low  flows;  i.e. 
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low  flows  in  the  winter  and  summer  may  belong  to  distinct  populations.  A  stronger 
possibility  is  that  the  available  low  flow  data  sets  do  exhibit  trend  and  non-random 
characteristics,  which  could  possibly  be  attributed  to  slow  cyclic  change  in  groundwater  levels 
or  to  climatic  trends.  Additiorul  testing  was  beyond  the  scope  of  the  current  investigations. 
However,  hirther  studies  are  recommended  since  these  results  may  call  into  question  the 
basic  assumptions  underlying  application  of  the  extreme  value  analysis  technique  for  analysis 
of  low  flow  characteristics. 

Stations  not  passing  the  statistical  screening  were  not  removed  from  the  analysis  as  there 
would  be  too  few  stations  remaining. 

3^.3      Summary  of  Data  Analysis  and  Low  Flow  Characteristics 

Application  of  various  non-parcunetric  tests  was  undertaken  for  the  available  data  base  for 
the  first  time  in  Ontario.^  The  test  results  indicate  that  the  available  data  base  of  extreme 
low  flov^^  may  exhibit  some  trend  and  dependence  with  some  possibility  of  non-random 
characteristics.  Previous  widespread  application  of  the  tests  utilized  have  not  been  found  in 
the  literature  for  low  flows.  The  average  length  of  record  for  the  stations  analysed  in  these 
regions  is  26  years.  It  is  possible  that  the  available  record  lengths  are  too  short  to  permit 
reasonable  application  and  interpretation  of  these  non-parametric  test  results.  Another 
possibility,  which  should  be  investigated  in  more  detail,  is  that  the  available  low  flow  data 
sets  do  exhibit  trend  and  non-random  characteristics.  The  latter  could  possibly  be  attributed 
to  slow  cyclic  duinges  in  groimdwater  levels,  seasonal  effects,  or  climatic  trends.  This  should 
be  investigated  in  more  detail  since  the  underlying  assumptions  for  application  of  the  extreme 
value  distribution  eind  subsequent  regression  analyses  are  called  into  question. 

The  Gumbel  Extreme  Value  Distribution  was  generally  found  to  adequately  fit  the  available 
low  flow  series  for  various  low  flow  durations.  However,  data  for  some  statioris  could  not 
be  fit  due  to  either  a  large  number  of  zero  flows,  numerical  computational  problems,  or  very 
large  low  flows  recorded  at  some  locations.  Additional  research  investigations  are  required 
to  develop  suitable  low  flow  analysis  techniques  for  such  locations. 


The  analysis  was  done  in  conjvinction  with  "Low  Flow  Characteristics  in  Ontario  Study", 
Cumming  Cockbum  Limited,  1988 
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Extreme  value  analyses  were  undertaken  on  an  annual  basis  for  125  stations.  A  total  of  73 
stations  were  retained  for  analysis  and  21  for  additional  testing  of  results. 

The  data  analyses  were  undertaken  for  both  regulated  and  unregulated  data  series. 
Therefore,  care  should  be  taken  in  comparison  and  interpretation  of  results,  notably  for  data 
series  which  may  include  the  effects  of  regulation. 

Figures  3.1  and  3.2  summarized  the  following  low  flow  characteristics: 

7  day  extreme  values  for  the  2,  5, 10  and  20  year  recurrence  intervals 

flows  which  were  equalled  or  exceeded  over  the  available  period  5, 50,  75,  95  and  99 

percent  of  the  time. 

Data  analysis  and  management  techniques  are  now  available  which  would  allow  efficiently 
updating  the  present  analyses  on  a  frequent  basis.  In  our  opinion,  the  low  flow  analyses 
should  be  updated  every  five  years  in  order  to  provide  reasonably  accurate  information  for 
investigation  requiring  low  flow  information. 

3.3         Physiographic  and  Meteoroloeic  Data 

3.3.1      Criteria 

Several  criteria  were  used  to  identify  physiographic  and  climatic  parameters  which  might  be 
suitable  for  use  in  regiorulizing  low  flow  characteristics.  These  are  discussed  as  follows: 

i)  Statistical  Significance 

When  undertaking  a  multivariate  analysis,  the  variables  chosen  must  make  a  contribution  to 
explaining  the  variance  of  the  low  flows.  The  significant  experience  in  undertaking  similar 
investigations  (discussed  in  Section  2.0)  was  used  to  identify  parameters  which  have  proven 
to  be  statistically  significant  in  predicting  low  flows. 
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ii)  Physical  Characteristics 


Wherever  possible,  variables  should  be  selected  based  on  hydrologic  significance.  That  is, 
the  parameters  should  have  some  physical  meaning  with  regard  to  estimates  of  low  flows. 

iii)  Reliability  of  Computation 

It  is  preferable  to  select  parameters  which  can  easily  be  computed  in  a  reliable  manner  by 
users  who  may  not  be  familiar  with  regression  procedures  or  the  details  of  the  statistical 
concepts.  Therefore,  from  a  practical  point  of  view,  it  was  desirable  to  make  the  parameter 
estimation  procedure  as  uncomplicated  as  possible  in  order  to  minimize  computation  errors 
when  applying  the  technique. 


3.3.2      Parameters 


The  parameters  selected  for  use  in  this  study  are  listed  as  follows: 


Hvdrometeorologic  Data 


Svmbol 


Index  of  mean  annual  precipitation  at  gauge  location  (mm) 
Index  of  mean  annual  snowfall  at  gauge  location  (cm) 
Index  of  mean  annual  runoff  at  gauge  location  (mm) 
Index  of  mean  annual  evaporation  at  gauge  location  (mm) 


MAP 
MAS 
MAR 
EVA 


Physiographic  Data 


Drainage  area  (km^) 

Index  of  area  controlled  by  lakes  and  swamps 

Length  of  main  channel  (km) 

10/85  channel  slope  (m/m) 

Base  Flow  Index  (dimensionless  ratio) 

Regulation  Index  (0  -  natural,  1  -  regulated) 

Groundwater  fluctuation  (m) 


DA 

ACLS 

LNTH 

SLP 

BH 

RN 

GW 
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Detailed  paraineter  definition  and  methodology  for  derivation  is  discussed  in  the  following 
sectioris: 

DAGcm^) 

The  watershed  drainage  area  was  obtained  from  records  published  by  the  Water  Survey  of 
Canada. 

ACLS  (%) 

An  index  representing  the  percentage  of  the  drainage  area  controlled  by  lakes  and  swamps 
(ACLS)  was  obtained  from  records  available  in  Regional  Flood  Frequency  Analysis  (Moin  and 
Shaw,  1986)  and  published  by  Environment  Canada. 

LNTHOan) 

The  length  of  the  main  channel  was  recorded  from  the  Regional  Flood  Frequency  Analysis 
(Moin  and  Shaw,  1986)  and  published  by  Environment  Canada.  The  lengths  for  some 
stations  were  not  previously  calculated.  For  these  stations  the  lengths  were  scaled  from 
review  of  Water  Survey  of  Canada  watershed  boundaries  and  rivers  on  1:50X)00  NTS  maps. 

SLOPE  (m/km) 

The  10/85  slope  (SLP)  recommended  by  the  U.S.  Geological  Survey  was  available  for  some 
stations  from  the  published  literature  Regional  Hood  Frequency  Analysis  (Moin  and  Shaw, 
1986).  The  10/85  slope  is  calculated  by  taking  tiie  elevation  difference  at  the  10%  and  85% 
points  along  the  channel  length.  The  slopes  for  some  stations  were  not  previously  calculated. 
For  these  stations  the  slopes  were  determined  by  creating  an  elevation  distance  table  from 
review  of  Water  Survey  of  Canada  watershed  boundaries  on  1:50,000  NTS  maps. 

Base  Flow  Index  CBFl)  (Dimensionless  Ratio) 
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the  drainage  basin.  BFI  is  defined  as: 

BFI  =  Total  Volume  of  Base  Flow 
Total  Volume  of  Runoff 

The  median  values  calculated  for  all  Ontario  gauging  stations  having  at  least  2  years  of 
continuous  daily  discharge  data  were  plotted  at  the  corresponding  drainage  basin  centroids 
using  1:600,000  scale  base  maps  for  Southern  Ontario  and  isoUnes  drawn.  The  centroids  were 
located  by  eye  after  delineating  the  draiiiage  area.  The  isoline  map  was  prepared  to  help 
provide  estimates  of  BFI  for  ungauged  basins  when  applying  the  regression  equations 
(Regional  Flood  Frequency  Analysis,  Moin  and  Shaw,  1986).  All  estimates  of  BFI  from  the 
isoline  maps  must  be  made  by  first  locating  the  basin  centroid  and  then  projecting  this  point 
to  the  closest  point  on  the  main  channel.  The  BFI  is  then  interpolated  from  the  isolines  at  this 
location  on  the  main  channel.  A  better  BFI  estimate  will  be  obtained  for  large  basins  and  in 
areas  where  the  isolines  are  very  dose  together,  if  an  average  value  of  BFI,  weighted  by  the 
area  between  isolines,  is  taken  over  the  entire  drainage  basin. 


RN 


Regulation  code  (i.e.  an  assessment  of  the  degree  of  regulation  a  stream  is  determined  to 
have): 


0  -  Natural,  non-regulated 

1  -  Regulated 

NOTE:  Heavily  regulated  stations  were  not  included  in  the  analysis  and  simple  correlation 
analysis  indicated  insignificant  levels  of  correlation  with  the  observed  low  flows  for 
the  stations  remeiining. 

GW(m) 


An  index  for  the  groundwater  table  fluctuation  was  developed  using  the  data  reported  in  the 
Enviroi\ment  Atks  of  Observation  WeUs  in  Ontario  published  by  the  Ontario  Ministry  of  the 
Environment,  1980.  Only  sparse  well  observation  information  is  available  but  nevertheless 
whatever  was  available  was  used.  The  maximum  fluctuation  of  the  water  table  for  a  well 
was  calculated  from  the  observed  maximum  and  minimum  water  table  depths. 
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MAP  (mm) 

An  index  of  mean  annual  precipitation  was  developed  with  reference  to  available  publications 
(MNR,  1984,  Environment  Canada,  1978  and  Flood  Frequency  Analysis,  Moin  and  Shaw, 
1986)).  The  index  of  mean  annual  precipitation  is  interpolated  for  each  watershed  used  in 
the  analysis  at  the  gauge  location. 

MAS  (mm) 

An  index  of  mean  annual  snowfall  was  obtained  from  available  information  published  by  the 
Ministry  of  Natural  Resources  (MNR,  1984)  and  Fisheries  and  Environment  Canada  ,  1978. 
An  index  map  was  used  to  derive  the  snowfall  index  for  each  hydrometric  gauge.  This  index 
represents  total  annual  snowfall  for  each  watershed.  The  index  was  determined  at  the  gauge 
location. 

MAR  (mm) 

The  index  of  mean  annual  runoff  is  expressed  as  a  depth  of  water  averaged  over  the  drainage 
basin  area.  The  isolines  of  runoff  were  obtained  from  existing  information  published  by 
Enviroiunent  Canada  and  adopted  by  the  Ministry  of  Natural  Resources  (Sangal  and  KaUio, 
1977,  MNR,  1984  and  Flood  Frequency  Analysis,  Moin  and  Shaw,  1988).  The  mean  aimual 
nmoff  index  was  derived  for  each  hydrometric  station  at  the  gauge  location. 

EVA  (mm) 

An  index  for  the  mean  annual  evaporation  was  obtained  from  the  publications  by  Fisheries 
and  Environment  Canada,  1978.  The  mean  annual  evaporation  index  was  derived  for  each 
hydrometric  station  at  the  gauge  location. 

3.3.3      Summary  of  Physiographic  and  Meteorologic  Data 

Tables  3.3(a)  and  3.3(b)  summarize  the  parameters  which  were  determined  for  each  station 
remaining  after  screening  for  the  Central  and  Southeastern  Regions,  respectively.  Table  3.4 
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3.3.3      Summary  of  Physiographic  and  Meteorologic  Data 

Tables  3.3(a)  and  3.3(b)  summarize  the  parameters  which  were  determined  for  each  station 
remaining  after  screerung  for  the  Central  and  Southeastern  Regions,  respectively.  Table  3.4 
summarizes  the  mean,  range  and  other  simple  statistics  of  the  data  base  for  the  Central, 
Southeastern  and  combined  regions. 


TABLE  3J  (a) 
DATA  BASE  FOR  THE  CENTRAL  REGION 


STATION 

GW 

RN 

YRS 

MAP 

MAS 

MAP 

SLP 

EVA 

DA 

BFI 

LNTH 

AHn^ 

702 

7O20 

n2EB008 

0.910 

1 

41 

990 

290 

.'?1(? 

0.020 

700 

1390 

0.79 

48.8 

100 

4.508 

1.459 

O2EB013 

0.910 

0 

46 

low 

220 

m 

4.010 

66,5 

593 

0.68 

38.0 

30 

4.738 

1.360 

02EC002 

0.910 

Q 

67 

89Q 

270 

454 

1.020 

740 

1520 

0.69 

94.3 

0 

1.690 

0.720 

02EC008 

0.460 

1 

10 

880 

170 

237 

3.302 

785 

274 

0.59 

30.0 

3 

0.267 

0.092 

02EC009 

0.300 

1 

22 

850 

160 

228 

1.314 

800 

181 

0.45 

29.0 

70 

0.293 

0.149 

02EC011 

0.940 

0 

16 

820 

230 

323 

0.420 

780 

282 

0.46 

55.3 

74 

0.250 

0.120 

02EC012 

0.910 

1 

13 

860 

240 

254 

1.640 

780 

324 

0.50 

32.5 

95 

0.360 

0.280 

02EC101 

0.790 

1 

12 

810 

200 

468 

4.610 

795 

24 

0.82 

32.0 

4 

0.250 

0.180 

02ED003 

0.370 

0 

34 

870 

260 

255 

1.900 

790 

1180 

0.55 

61.8 

51 

2.260 

1.450 

02ED0O5 

0.430 

1 

19 

880 

270 

424 

2.230 

785 

295 

0.61 

50.5 

35 

0.970 

0.570 

02ED007 

0.910 

0 

17 

1000 

300 

408 

1.990 

740 

177 

0.63 

28.0 

0 

0.960 

0.780 

02ED009 

0.150 

Q 

11 

890 

270 

295 

1.620 

76^ 

95 

0.33 

15.5 

23 

0.030 

0.010 

02ED011 

0.910 

1 

]4 

980 

220 

438 

2.177 

72.5 

168 

0.65 

14.0 

5 

0.463 

0.217 

02ED100 

0.270 

1 

14 

820 

230 

209 

2.750 

800 

86 

0.57 

17.8 

0 

0.150 

0.060 

02ED103 

0.380 

Q 

14 

850 

260 

380 

4.820 

790 

195 

0.71 

37.5 

0 

0.880 

0.700 

02HB(X)1 

3.050 

1 

58 

810 

120 

271 

1.800 

810 

205 

0.45 

15.0 

0 

0.603 

0.317 

O2HB0O5 

0.550 

1 

24 

780 

130 

394 

3.360 

825 

96 

0.53 

16.5 

0 

0.230 

0.060 

02HB008 

0.270 

1 

20 

790 

160 

324 

4.290 

815 

127 

0.57 

21.0 

0 

0.360 

0.230 

02HB011 

0.370 

1 

19 

790 

140 

361 

2.870 

830 

235 

0.64 

30.0 

3 

0.450 

0.330 

02HB013 

3.050 

1 

20 

870 

140 

325 

5.927 

800 

62 

0.65 

9.0 

40 

0.211 

0.144 

02HC(X)5 

0.080 

1 

35 

790 

140 

350 

3.810 

815 

88 

0.40 

18.0 

0 

0.160 

0.058 

02HC006 

0.080 

1 

29 

790 

150 

349 

7.470 

810 

249 

0.54 

29.5 

38 

0.740 

0.370 

02HC012 

0.430 

1 

25 

780 

200 

271 

5.350 

810 

169 

0.60 

24.6 

24 

0.480 

0.310 

02HC013 

0.060 

0 

17 

780 

140 

3^S 

6.000 

810 

88 

0.43 

15.2 

0 

0.230 

0.080 

02HC017 

1.219 

0 

19 

790 

160 

285 

3.378 

820 

63 

0.26 

20.3 

0 

0.046 

0.003 

02HC018 

0.183 

0 

24 

800 

140 

258 

7.870 

810 

106 

0.46 

21.5 

2 

0.104 

0.042 

02HC019 

0.610 

Q 

24 

790 

140 

402 

9.420 

810 

94 

0.58 

17.5 

7 

0.461 

0.327 

02HC022 

0.457 

1 

26 

800 

150 

239 

2.910 

810 

186 

0.37 

28.1 

13 

0.138 

0.069 

02HC023 

0.366 

0 

25 

790 

180 

242 

4.438 

820 

62 

0.57 

11.8 

5 

0.183 

0.125 

02HC024 

0.061 

0 

25 

780 

140 

380 

3.563 

820 

316 

0.49 

37.5 

3 

1.387 

1.162 

02HC025 

0.366 

0 

25 

780 

170 

254 

1.571 

820 

303 

0.59 

45.0 

5 

0.820 

0.591 

02HC028 

0.610 

1 

23 

780 

140 

314 

3.302 

810 

78 

0.34 

17.5 

0 

0.089 

0.045 

02HC029 

0.488 

1 

23 

780 

140 

347 

4.212 

820 

130 

0.43 

20.0 

0 

0.418 

0257 

02HC030 

0.610 

1 

21 

780 

140 

323 

4.410 

800 

204 

0.24 

38.0 

0 

0.228 

0.142 

02HC031 

0.488 

1 

18 

790 

180 

209 

4.168 

820 

148 

0.23 

24.0 

0 

0.009 

0.000 

02HC032 

0.549 

1 

21 

790 

180 

189 

1.841 

810 

95 

0.44 

21.0 

19 

0.071 

0.052 

02HC033 

0.488 

1 

21 

780 

140 

M2 

3.267 

820 

71 

0.25 

29.5 

0 

0.096 

0.045 

02HC034 

0.488 

0 

18 

780 

140 

202 

6.190 

820 

194 

0.17 

22.0 

0 

0.006 

0.000 

02HDOO3 

0.488 

1 

28 

800 

140 

527 

3.178 

790 

67 

0.70 

15.5 

0 

0.484 

0.299 

02HDOO6 

0.488 

1 

2« 

810 

140 

483 

5.680 

790 

83 

0.59 

22.0 

7 

0.503 

0.353 

02HD009 

0.488 

0 

22 

810 

140 

359 

1.153 

790 

83 

0.62 

15.0 

0 

0.341 

0.228 

02HD010 

0.488 

0 

22 

810 

140 

408 

8.126 

790 

65 

0.62 

17.5 

0 

0.316 

0.195 

02HD012 

0.488 

0 

11 

810 

160 

449 

9.434 

790 

232 

0.66 

21.0 

0 

1.388 

1.036 

02HH002 

0.549 

1 

15 

900 

180 

416 

2.540 

750 

326 

0.65 

30.0 

100 

0.702 

0.177 

02HJ001 

0.488 

0 

25 

790 

180 

322 

2.714 

780 

110 

0.40 

18.5 

89 

0.076 

0.025 

02HJ003 

0.366 

1 

20 

790 

170 

320 

1.065 

780 

282 

0.62 

28.5 

77 

0.156 

0.023 

02HKOO'; 

0.427 

1 

19 

910 

180 

46? 

0.920 

770 

456 

0.76 

35.5 

87 

0.803 

0.226 

02HK006 

0.427 

1 

14 

900 

170 

424 

2.743 

770 

541 

0.62 

60.0 

30 

0.157 

0.021 

TABLE  3  J  (b) 
DATA  BASE  FOR  THE  SOUTHEASTERN  REGION 


STATION 

GW 

RN 

YR? 

MAP 

MAS 

MAR 

SLP 

EVA 

DA 

BFl 

LNTH 

ACLS 

76i 

■  "i^-iT 

02HE002 

2.13 

13 

850 

170 

417 

0.880 

790 

114 

0.32 

15.0 

95 

0.002 

0.000 

02HK005 

0.91 

14 

910 

180 

462 

0.920 

765 

456 

0.76 

35.5 

87 

0.803 

0.226 

02HK006 

1.16 

13 

920 

170 

353 

3.441 

760 

541 

0.62 

31.0 

5 

0.157 

0.021 

02HL001 

2.74 

67 

810 

170 

361 

0.960 

780 

2620 

0.71 

87.0 

56 

1.477 

0.582 

02HLO03 

2.74 

18 

850 

170 

384 

2.200 

770 

401 

0.60 

66.5 

69 

0.450 

0.061 

02HM002 

4.27 

30 

880 

180 

375 

7.144 

770 

189 

0.65 

20.0 

100 

0.322 

0.039 

02HM003 

4.27 

24 

850 

180 

366 

1.530 

780 

891 

0.77 

112.4 

79 

0.145 

0.034 

02HM005 

4.27 

0 

10 

900 

180 

540 

0.960 

775 

155 

0.49 

27.0 

35 

0.025 

0.000 

02HM006 

4.27 

13 

900 

18Q 

462 

0.800 

775 

150 

0.63 

40.0 

84 

0.095 

0.039 

02HM007 

4.27 

10 

880 

180 

403 

0.610 

770 

694 

0.78 

50.5 

97 

1.159 

0.355 

02KA003 

0.30 

0 

20 

830 

190 

260 

12.190 

635 

7 

0.50 

5.0 

50 

0.001 

0.000 

02KB  001 

0.30 

67 

800 

190 

362 

1.360 

660 

4120 

0.88 

107.6 

59 

11.292 

6.560 

02KC009 

0.30 

61 

800 

190 

355 

3.217 

715 

2380 

0.68 

90.0 

50 

4.762 

2.830 

02KF011 

0.37 

0 

12 

810 

200 

358 

0.280 

720 

269 

0.33 

35.0 

2 

0.059 

0.010 

02KF012 

0.61 

11 

810 

200 

375 

1.720 

720 

203 

0.61 

37.5 

61 

0278 

0.181 

02KF013 

0.85 

11 

800 

200 

387 

1.280 

735 

280 

0.67 

32.5 

45 

0.150 

0.021 

02KF014 

0.55 

12 

800 

180 

358 

1.725 

740 

277 

0.59 

53.0 

95 

0.109 

0.000 

02LA004 

0.61 

34 

810 

200 

m 

0.940 

735 

3830 

0.52 

92.5 

47 

5.705 

3.118 

02LA006 

0.61 

13 

820 

200 

412 

0.140 

740 

409 

0.41 

63.0 

63 

0,076 

0,008 

02LA007 

0.61 

0 

14 

850 

190 

380 

0.610 

740 

559 

0.38 

60.5 

17 

0.086 

0.012 

02LB006 

1.83 

0 

35 

820 

190 

420 

0.650 

740 

433 

0.30 

36.3 

0 

0.195 

0.089 

02LB0O8 

0.30 

0 

25 

830 

190 

478 

0.540 

730 

440 

0.30 

46.2 

0 

0.161 

0.105 

02LB017 

0.91 

0 

9 

800 

200 

395 

1.016 

750 

69 

0.30 

15.0 

0 

0.529 

0.097 

02LB022 

2.44 

0 

9 

800 

200 

396 

2.540 

745 

152 

0.30 

34.0 

0 

0.027 

0.001 

02MCXX)1 

0.30 

0 

26 

800 

200 

422 

0.293 

745 

404 

0.37 

26.0 

0 

0.044 

0.004 
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4.0  ANALYSES  OF  EXISTING  TECHNIQUES 

4.1  General 

The  previous  results  (Southwestern  and  West  Central  Regions,  CCL,  1990)  indicated  that: 

•  estimates  of  low  flow  characteristics  could  not  be  improved  by  subdividing  the  area 
into  sub-regions  (perhaps  due  to  the  relatively  small  amount  of  available  data),  and 

•  in  regard  to  estimating  accuracy,  the  methods  ranked  as  follows: 

i)  Regression 

ii)  Mapped  IsoUnes 

iii)  Index  Relationships 

iv)  Nearby  gauge  proration 

To  assess  the  appUcabUity  of  the  previous  results  extrapolated  to  the  present  study  area,  the 
previously  developed  techniques  with  the  previously  determined  coefficients  and  constants 
were  utilized  to  predict  low  flow  characteristics  for  stations  located  in  the  Central  and 
Southeastern  Regions. 

4.2  Summary  of  Analysis 

This  analysis  is  summarized  in  Table  4.1  for  the  Central,  Southeastern  and  the  combined 
regions.  (A  detailed  analysis  is  given  in  Appendix  E.)  The  results  indicate  that  the  ranking 
of  the  technique  remains  almost  the  same  as  the  previous  test  stations  indicated 
[(i.e.  i)  Regression,  ii)  Index  Relationships,  iii)  Mapped  IsoUnes,  iv)  Unit  Area  Proration 
techniques.] 

The  Nash-Sutdiffe  R^  values  were  computed  and  indicated  similar  results  for  application  of 
the  regression  equation  techniques  compared  to  the  previous  study  (i.e.  R^  for  estimation  of 
7Q20  ranging  from  0.29,  0.75  and  0.70  for  the  Central,  Southeastern  and  combined  region 
analysis).  In  addition,  the  Index  method  and  Isoline  method  are  reversed  in  ranking.  This 
could  be  due  to  a  different  relationship  between  BFI  and  unit  area  low  flows  for  these 


TABLE  4.1 

SUMMARY  OF  TESTING  PREVIOUS  LOW  FLOW  ESTIMATION  TECHNIQUES 

NASH-SUTCLIFFE  R' 

FOR  REGIONAL  ANALYSIS 


Estimating  Method  and 
Target  Characteristic 

Central 

Southeastern 

Combined 

Regression  Equation 

70: 

0.24 

0.88 

.75 

7Q20 

0.29 

0.75 

.70 

Index  Method 

70: 

0.22 

0.59 

.53 

7Q20 

0.04 

0.42 

.39 

Isoline  Method 

7Q20 

-0.09 

-0.25 

-020 

Station  Proration 

7Q20 

-0.94 

-1.41 

-129 
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regions.  It  is  interesting  to  note  that  in  all  regions  the  unrefined  techniques  still  produce 
better  results  than  nearby  gauge  proration  techniques  historically  used. 

Table  4.2  summarizes  the  average  7Q20  low  flow  rates  and  standard  deviations  for  the 
observed  (i.e.  statistically  determined  flows)  and  estimation  techniques  for  each  region. 

The  Index  method  produced  lower  average  low  flows  for  all  three  regional  groups.  This 
indicates  that  drainage  areas  in  these  regions  produce  more  low  flow  per  unit  area  compared 
to  the  values  in  Southwestern  and  West  Central  Ontario  from  which  the  index  curves  were 
developed. 

The  nearby  gauge  proration  technique  and  isoline  technique  both  produced  average  values 
greater  than  the  observed  values  for  these  regions.  This  preliminary  isoline  analysis  was 
based  on  conversion  of  BFI  isolines  to  unit  area  low  flows  (i.e.  the  comparison  was  developed 
from  data  from  the  Southwestern  and  West  Central  regions).  This  would  appear  to  indicate 
that  BFI  vedues  in  the  Central  and  Southeastern  region  produce  higher  low  flows  than  similar 
BFI  values  in  the  Southeastern  and  West  Central  regions.  Conversely,  the  low  flow /BFI 
relatioiiship  may  be  being  affected  by  Basin  Storage  BFI  interrelationships  (see  Section  2.1). 

The  regression  technique  estimated  lower  average  low  flows  in  the  Central  region  and  the 
two  regions  combined.  However,  higher  average  low  flows  were  predicted  for  the 
Southeastern  Region.  This  suggests  that  the  Central  and  Southeastern  regions  are  likely  to 
be  distinct  from  Southwestern,  West  Central  Ontario  and  distinct  from  each  other. 

In  most  cases,  estimation  techniques  for  the  low  flows  with  smaller  recurrence  intervals  (i.e. 
7Q2)  resulted  in  higher  Nash-Sutcliffe  R'  values  than  did  the  low  flows  with  less  frequent 
recurrence  intervals  (i.e.  7(22o)-  Therefore,  it  might  be  useful  to  consider  fliat  further 
investigation  should  examine  using  ZQj  as  the  predictor  and  then  some  relatior\ship  to 
produce  7Q20  values. 

Summary 

•        Nash-Sutcliffe  R^  values  are  quite  reasonable  for  the  regression  equations  for 
Southeastern  and  the  combined  regions. 


TABLE  iJZ 

CHARACTERISTICS  OF  OBSERVED 

AND  SIMULATED  7Qa  VALUES 

FROM  PREVIOUS  TECHNIQUES 


Techniques 

Central 

Southeastern 

Combined 

Average 
7Q^ 

Standard 

Deviation 

7Q« 

Average 

7Q« 

Standard 

Deviation 

7Q2C 

Average 
7Q» 

Standard 

Deviation 

7Qa, 

Observed 

.31 

.39 

.64 

154 

.41 

.94 

Regression 

.22 

.20 

1.0 

1.9 

.47 

1.14 

Index 

.11 

.10 

0.4 

0.47 

.19 

.30 

isoline 

.38 

.56 

13 

2.8 

.72 

1.73 

Unit 

.36 

.53 

0.9 

2.4 

.53 

1.46 
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None  of  the  estimation  techniques  produced  favourable  results  in  the  Central  Region. 

Low  flows  (i.e.  7Q20  are  lower  for  methods  based  only  on  drainage  area  and  higher 

for  methods  based  on  BFI  values.   Hence,  coefficients  for  these  parameters  as  they 

relate  to  low  flows  in  the  Central  and  Southeastern  regions  are  probably  significantly 

different  than  from  the  Southwestern  West  Central  regions  whence  they  were 

developed. 

The  large  difference  in  Nash  Sutcliffe  R^  values  for  the  Central  and  Southeastern 

regior\s  indicate  they  may  be  distinct  with  respect  to  low  flow  characteristics 

The  isoline  method  based  upon  BFI  values  did  not  appear  to  provide  a  good 

estimation  technique 

The  Nash- Sutcliffe  R^  results  indicate  that  all  techniques  are  better  predictions  than 

use  of  the  gauge  unit  area  proration  technique  which  is  in  common  use. 
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5.0  REGIONALIZATION  OF  LOW  FLOW  CHARACTERISTICS 

5.1  General 

Four  alternative  methods  of  providing  regional  estimates  of  low  flow  characteristics  were 
considered.  They  are;  Multiple  Linear  Regression  (Section  5.2),  Index  Method  (Section  5.3), 
Mapped  Isolines  of  imit  low  flows  (Section  5.4),  and  proration  from  nearby  gauges  (unit  area 
technique)  (Section  5.5).  Emphasis  was  placed  on  determining  7Q20  and  ZQ,  as  these  flows 
were  identified  to  be  the  key  low  flow  statistics  required  by  the  Ministry  of  the  Environment. 

5.2  Reeression 
5.2.1     General 

Multiple  linear  regression  equations  take  the  following  general  form: 

Y  =  ao  +  a,Yi+  ajYj  + +  a„Y„  (5.1) 

where  Y  =    the  dependent  variable  (e.g.  ZQ^q) 

Yj  Yj.-.Yn  =      the  independent  variables  (e.g.  physiographic  and  meteorologic 

watershed  characteristics) 
ao,ai...a^  =       regression  coefficients 

In  order  to  obtain  a  more  accurate  regression  equation,  it  is  sometimes  necessary  to  transform 
the  data  by  taking  logarithms,  square  roots,  cubes,  etc.  The  transformations  considered  for 
this  investigation  are  discussed  in  Section  5.2.3. 

A  number  of  regression  procedures  are  available,  including  development  of  all  possible 
equations,  forward  selection,  backward  selection,  stagewise  regression  and  stepwise 
regression.  The  stepwise  regression  procedure  is  generally  recommended  for  use  in  practical 
applicatior^s  (Draper  and  Smith,  1981)  and  was,  therefore,  adopted  for  the  purposes  of  this 
investigation. 
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5^^      Methodology 

The  regression  equations  were  develop>ed  utilizing  the  stepwise  multiple  Unear  regression 
procedure  available  in  the  Statistical  Package  for  the  Social  Sciences  (SPSS)  (Nie,  1975).  More 
specifically,  the  regression  subprogram  has  been  used  and  has  the  following  special  features: 

1)  Out  of  various  procedures  for  selecting  variables,  including  forward  selection, 
backward  elimination  and  the  stepwise  selection,  the  last  one  which  is  really  a 
combination  of  backward  and  forward  procedures,  was  adopted  for  this  investigation 
as  it  is  the  most  commonly  used  procedure  for  this  type  of  study. 

2)  Variable  selection  -  all  independent  variables  can  be  stored  and  then  only  those 
variables  desired  for  a  particular  analysis  called  up  and  used  according  to  the  desired 
form  of  the  equation 

3)  Combination  of  variables  -  variable  transformation  and  new  variables  may  be 
computed  from  existing  variables 

4)  Transformations  -  the  variables  may  be  transformed  (e.g.  square  root,  logarithmic, 
squared,  etc.)  in  order  to  more  nearly  linearize  the  relationships 

5)  Calculation  of  statistics  -  the  SPSS  regression  package  allows  calculation  of  numerous 
regression  coefficients,  statistics  and  residual  statistics  (difference  between  observed 
and  calculated  values).  Also  possible  are  scatter  plots  of  residuals  and  statistical  tests 
for  residual  analysis,  etc. 

The  regression  constant  and  regression  coefficients  are  determined  in  order  to  minimize  the 
sum  of  the  square  residuals.  The  residuals  are  the  difference  between  the  observed, 
dependent  variable  and  the  prediction  by  the  regression  equation.  The  SPSS  program 
automatically  includes  those  independent  variables  which  meet  the  95%  confidence  level 
based  on  the  computed  values  of  the  F  statistic.  Only  those  variables  which  meet  the 
specified  level  at  any  given  step  are  retained  in  the  regression  equation  and  all  those  variables 
which  fall  below  the  specified  level  are  deleted  from  the  regression  equation.  By  specifying 
an  appropriate  subcommand,  a  variable  can  be  forced  in  or  out  of  the  regression  equation. 

The  entry  of  the  variables  can  also  be  controlled  by  specified  the  maximum  number  of  steps 
and  redefining  the  tolerance  criteria.  Variables  must  pass  both  tolerance  and  minimum 
tolerance  tests  in  order  to  enter  and  remain  in  a  regression  equation.    Tolerance  is  the 


5-3 

propwrtion  of  the  variance  of  a  variable  in  the  equation  that  is  not  accounted  for  by  other 
independent  variables  in  the  equation.  The  minimum  tolerance  of  a  variable  not  in  the 
equation  is  the  smallest  tolerance  any  variable  already  in  the  equation  would  have  if  the 
variable  being  considered  v^ere  included  in  the  analysis.  If  a  variable  passes  the  tolerance 
criteria,  it  is  further  tested  according  to  the  method  in  effect. 

5.2.3      Transformed  and  Derived  Parameters 

The  transformations  used  for  this  analysis  were  logio,  square  root,  square  and  cube  for  the 
meteorologic  and  physiographic  parameters.  The  derived  parameters  used  in  this  analysis 
are  defined  below: 

1.    S>U9)e  Factor  1,  SFI  =      ^^ 


2.    Shape  Factor  2,  SF2  = 


LNTH^ 

DA 
LNTH 


3.    Shape  Factor  3,  SF3  =  DA*SLP* 


LNTH 


2 

4.  Infiltration  or  Recharge  or  Deep  Percolation,  INF  =  MAP -MAR-EVA 

5.  PI  =  {LNTHI  yfsil^'* 

6.  P2  ^  (BFII  (LNTH/  y/SLP)"^ 

The  derived  p«u-ameters  SFI,  SF2,  SF3  and  INF  were  newly  defined  specifically  for  this 
investigation. 

5.2.4      Simple  Correlation  of  Parameters 

Simple  correlations  between  independent  parameters  and  between  independent  and 
dependent  parameters  were  examined  to  screen  parameters  for  input  to  subsequent 
regression  analyses.  The  matrices  of  simple  correlation  coefficients  for  the  Central, 
Southeastern  and  combined  regions  are  presented  in  Appendix  F.  A  significant  correlation 
coefficient  for  a  95%  level  of  significance  has  a  magnitude  equal  or  greater  than  0.282,  0.388 
and  0.226  for  each  of  the  Central,  Southeastern  and  combined  regions  respectively.  It  is 
evident  from  the  tables  that  a  few  parameters  exhibit  significant  intercorrelation. 
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In  general  it  was  found  that  the  Regulation  Code  (RN)  is  not  significantly  correlated  with  the 
low  flows  except  for  the  Central  region  for  7Q20  and  the  magnitude  of  the  correlation 
coefficient  there  also  is  small. 

It  was  also  found  that  the  independent  parameters  which  are  generally  most  highly  correlated 
to  the  dependent  parameters  are  DA,  EVA  and  BFI.  (Some  other  parameters  were  found  to 
be  inter-correlated  with  these.)  Some  additional  data  transforms  were  undertaken  (see 
Appendix  A)  which  indicated  that  the  cube  of  the  drainage  area  remained  highly  correlated 
with  the  dependent  paramieters. 

5.2.5      Regression  Equation  Development 

A  large  number  (several  hundrtxi)  preliminary  regression  equations  were  developed  in  order 
to  predict  the  7Q2  and  ZQjo  low  flows  as  a  function  of  basin  physiographic  and 
hydrometeorologic  parameters  (the  procedure  is  described  in  Appendix  A). 

The  final  form  of  the  equations  developed  are  summarised  in  Tables  5.1  and  5.2.  The  final 
form  of  equation  selected  takes  into  account  the  prediction  robustness  assessed  during  testing 
on  the  reserved  stations  (see  Appendix  A). 

A  sioitable  level  of  accuracy  for  the  Southeastern  region  could  not  be  achieved  without 
inclusion  of  the  cube  of  drainage  area  as  an  independent  parameter.  This  leads  to  a  small 
regression  coefficient  (aj),  and  while  this  is  undesirable  from  a  practical  applications 
viewpoint,  such  coefficients  were  found  to  be  statistically  valid.  Nevertheless,  in  our  opinion, 
additional  research  is  required  to  confirm  that  the  Southeastern  Region  is  imusual  with 
resp)ect  to  low  flows  and  to  investigate  possible  alternative  forms  of  the  prediction  equation. 

In  addition,  through  analysis  of  interim  results,  it  was  noted  that  the  Central  and 
Southeastern  regions  were  distinct  regions  with  regard  to  their  low  flow  characteristic 
interrelationships  with  the  physiographic  and  climatic  parameters.  This  confirms  that  it  is 
not  appropriate  to  analyse  them  as  combined  regions  for  the  regression  technique. 


TABLE  5.1 
SUMMARY  OF  REGRESSION  ANALYSIS  (CENTRAL  REGION) 


Y=aO+al*DA+a2*Bn 


Dependent 
Parameter 

Independent  Parameters 

aO 

al 

a2 

N 

SE 

R^2 

Central  R 

egion 

7O20 

-0.2134 

6.6184E-04 

0.7022 

48 

0.29 

0.53 

702 

-0.7216 

1.8060E-03 

1.7386 

48 

0.73 

0.55 

302 

-0.5398 

1.6260E-03 

1.2856 

48 

0.45 

0.70 

3O20 

-0.1841 

5.8893E-04 

0.6295 

48 

0.27 

0.51 

3O50 

-0.1331 

4.5199E-04 

0.5160 

48 

0.25 

0.42 

30O2 

-0.7119 

2.2380E-03 

1.6806 

48 

0.63 

0.70 

30O20 

-0.3275 

9.7749E-04 

0.9305 

48 

0.36 

0.59 

30O50 

-0.2839 

8.7086E-04 

0.8045 

48 

0.33 

0.57 

1 

TABLE  52 
SUMMARY  OF  REGRESSION  ANALYSIS  (SOUTHEASTERN  REGION) 


Y=aO+al*DA^3+a2*Bn 


Dependent 
Parameter 

Independent  Parameters 

aO 

al 

a2 

N 

SE 

R'^2 

Southeastern  Region 

7O20 

-0.5084 

7.6323E-11 

1.1460 

25 

0.64 

0.88 

702 

-0.9018 

1.3049E-10 

2.2728 

25 

0.98 

0.91 

302 

-1.0351 

1.2409E-10 

2.3828 

25 

0.93 

0.91 

3O20 

-0.6133 

7.0980E-11 

1.2527 

25 

0.64 

0.89 

3O50 

-0.6226 

6.5153E-11 

1.2372 

25 

0.64 

0.85 

30O2 

-1.0195 

1.4637E-10 

2.6144 

25 

1.03 

0.92 

30O20 

-0.5196 

8.5495E-11 

1.3062 

25 

0.74 

0.88 

30O50 

-0.4643 

7.9836E-11 

1.1773 

25 

0.70 

0.87 
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5.3        Index  Method 

The  regression  analysis  confirmed  the  conclusions  drawn  from  the  literature  survey  that 
draiiuge  area  (DA)  is  a  good  predictor  of  low  flows.  Therefore,  a  simple  method  using  DA 
alone  to  estimate  low  flows  was  investigated.  To  meet  this  objective,  graphs  of  VQj  and  7Q20 
as  a  function  of  DA  are  plotted  on  Figure  5.1  for  the  Central,  Southeastern  and  the  combined 
region.  A  best  fit  analysis  was  performed  and  results  of  the  regression  are  summarized  in 
Table  5.3. 

The  ratios  of  7Qy  to  various  n  day  flows  with  y  year  recurrence  were  determined  and  plotted 
on  Figure  5.2.  Finally,  the  ratio  of  7Qy  to  7Q2  for  all  these  regions  was  determined  as  shown 
on  Figure  5.3. 

To  use  this  method  knowing  the  drainage  area,  the  7Q2  flow  is  estimated  from  Figure  5.1  (or 
equation  in  Table  5.4)  for  the  appropriate  region  in  which  the  ungauged  watershed  is  located. 
The  7Qy  flow  can  then  be  estimated  using  Figure  5.2  for  the  required  recurrence  interval  (y). 
Low  flows  for  other  n  day  durations  can  then  be  estimated  using  Figure  53. 

The  graphs  of  n-day  low  flows  to  7  day  low  flows  as  well  as  the  graph  of  7  day  y  year 
recurrence  ratio  to  7Q2  low  flows  give  some  insight  into  the  interrelationships  between  the 
various  extreme  value  low  flows.  For  example,  the  ratio  of  7  day  average  low  flows  to  30 
day  low  flows  tends  to  decrease  slightly  v^th  increased  recurrence  interval  while  the  ratio 
of  7  day  average  low  flows  to  3  day  average  low  flows  increases  with  increased  recurrence 
interval.  As  expected,  the  30  day  low  flows  are  greater  than  7  day  low  flows  (i.e.  ratio  ranges 
from  1.1  to  1.3  times,  1.5  to  1.9  times,  and  1.3  to  1.4  times  the  7  day  low  flow  values  for  the 
Centred,  Southeastern  and  combined  regions). 

The  development  and  use  of  the  Index  Method  indicates  that  it  may  be  possible  to  utilize 
short  periods  of  record  to  estimate  7Q2  flows  and  then  use  the  other  graphs  developed  here 
to  estimate  low  flows  for  other  recurrence  intervals  (i.e.  a  few  years  of  data  may  provide  a 
reasonable  estimate  of  7Q2  which  can  then  be  extended  by  use  of  the  graphs). 


TABLE  53 

SUMMARY  OF  BEST  FFT  ANALYSIS 

FOR  7Q,  AND  70^,  AS  A  FUNCnON  OF  DRAINAGE  AREA 

Y  =  Ao  +  A,  •  DA 


Dependent 
Parameter 

Independent  Parameters 

ao 

a, 

N 

SE 

R' 

Central  Region 

7Q» 

0209 

5.8910-' 

48 

035 

0.46 

7Q2 

0.383 

1.61*1 0-' 

48 

1.11 

038 

Southeastern  Region 

7Q» 

-1.008 

1.46*10^ 

25 

1.371 

0.7 

7Q2 

-1.60 

25110^ 

25 

2.182 

0.73 

Combined  Central  and  Southeastern  Regions 

7Q« 

-0.198 

1.180*10^ 

73 

0.98 

0.71 

7Q: 

0.118 

2.0510^ 

73 

1.72 

0.71 
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5.4        Mapped  Isoline  Method 

Previous  studies  have  used  interrelationships  between  mapped  isolines  of  BFI  which  had 
been  converted  to  relate  to  unit  area  low  flow  characteristics  (Gumming  Cockbum  Limited, 
1989).  Further  to  the  previous  analysis,  isolines  of  urut  area  7Q2  low  flows  were  drawn  for 
the  Central  and  Southeastern  combined  regions  (see  Figures  5.3  a,  b,  c  and  d)  (in  pocket) 
using  digital  terrain  modelling  procedures. 

First,  a  data  set  of  station  locatior^s  (x  and  y)  and  unit  area  low  flow  characteristics  (i.e.  2  = 
7Q2)  were  compiled.  Then  the  data  set  was  used  by  the  GWN  digital  terrain  modelling 
package  to  create  a  trianguiation  matrix  over  the  region  on  which  to  interpolate  the  location 
of  "even  values"  isoline  intersection  points.  The  isolines  are  then  created  based  on  the  above 
interpolation. 

The  density  of  isolines  is  a  reflection  of  station  density.  The  higher  density  in  the  Central 
region  results  in  closely  packed  isolines  while  the  Southeastern  region  isolines  are  more 
spread  out  due  to  lower  station  density. 

Several  patterns  seem  to  become  evident  using  this  procedure.  For  example,  there  appears 
to  be  some  kind  of  lake  effect  which  influences  islines  near  Lake  Ontario  and  some  influence 
of  the  Caixadian  Shield  on  isolines  in  proximity  to  the  Shield. 

In  addition,  a  map  of  imit  area  7Q20  isolines  was  produced  for  the  combined  regions.  As 
expected,  the  map  has  the  same  general  appearance  as  the  individual  region  7Q2  isoline  maps. 
However,  in  general,  there  are  fewer  stations  and  more  stations  had  ZQ^o  =  0  in  the  SE  region. 
Therefore,  7Q2  isolines  were  cor\sidered  more  appropriate  and,  for  this  study,  further 
investigations  used  oiUy  the  uiut  area  7Q2  mapped  isoUnes.  Future  investigations  should 
consider  using  the  7  day  mean  low  flow  as  an  index  value  and  calculating  this  value  for 
stations  with  short  term  records  (e.g.  5  years)  in  order  to  refine  the  7Q2  isoline  maps  in 
Southeastern  Ontario. 

Utilizing  the  combined  map  provides  isolines  based  on  a  greater  number  of  data  points  which 
provide  better  overall  results. 
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5.5  Station  Proration 

In  the  past  it  has  been  common  practice  to  prorate  unit  flows  from  nearby  gauged 
watersheds.  This  is  generally  done  by  experienced  hydrologists  who  have  a  good 
understanding  of  local  stream  characteristics  and  other  factors  within  the  region  (i.e. 
diversions  and  regulation,  etc.). 

For  assessment  of  this  method,  the  7Q20  low  flow  characteristics  were  determined  for  statioris 
in  the  region  and  svunmaiized  in  Figures  3.1  and  3.2  (in  pocket). 

To  estimate  low  flows  at  a  location  between  two  or  more  sites  located  on  Figures  3.1  and  3.2, 
the  reciprocal  distance  p)ercentage  is  used  in  the  proration  (i.e.  the  distance  to  nearby  gauges 
would  be  estimated  and  unit  values  are  then  weighted  by  the  reciprocal  of  the  distance  as 
a  p)ercentage  and  then  averaged). 

5.6  Summary 

Of  the  above  four  low  flow  estimating  methods  described  in  the  previous  sections.  Mapped 
Isolines  and  the  Index  Method  are  the  easiest  to  use.  The  Station  Proration  Method  has  been 
included  due  to  its  widespread  use  and  to  evaluate  whether  the  alternate  techniques  provide 
improved  estimation  of  low  flow  characteristics.  The  relative  prediction  accuracy  of  all  these 
methods  is  evaluated  and  discussed  in  Section  6.0. 
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6.0  TESTING  PREDICTION  METHODS 

6.1  General 

The  relative  prediction  accuracy  of  the  foiir  methods  develojjed  in  Section  5.0  was  examined 
using  12,  9  and  21  stations  for  the  Central,  Southeastern  and  the  combined  regions, 
respjectively.  These  stations  were  not  included  in  the  development  of  the  regional  low  flow 
estimation  techniques. 

6.2  History  of  Test  Stations 

The  test  stations  all  have  at  least  10  years  of  record,  and  are  active  within  the  last  5  years. 
Stations  were  selected  spatially  to  cover  the  complete  region  and  with  a  range  of  similar 
characteristics  to  the  remaining  data  base.  The  data  base  for  these  statioi\s  is  presented  in 
Tables  6.1  and  62  for  the  Central  and  Southeastern  regions,  respectively. 

Further  investigation  indicated  that  two  stations,  02KC014  and  02LB013  are  heavily  regulated 
such  that  low  flows  may  be  affected.  For  the  Southeastern  and  hence,  the  combined  Central 
and  Southeastern  regions,  the  analysis  was  conducted  by  including  stations  02KC014  (Indian 
River  near  Pembroke)  and  02LB013  (South  Nation  River  at  Casselman)  in  the  test  stations  and 
also  by  excluding  these  two  stations  from  the  test  stations. 

6.3  Goodness  of  Fit 

To  test  the  goodness  of  fit,  Cumming  Cockbum  Limited  (1990)  used  the  Nash-Sutdiffe  (1970) 
model  efficiency  (N.S.R.^).  The  N.S.R.^  is  calculated  by  the  following  relationship: 

Nsji'  - 1  -  IMsM  (1) 

where  Q^  and  Q,  are  the  observed  and  simulated  discharges,  and  Q„  is  the  mean  of  the 
observed  discharges. 
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The  "observed"  low  flows  (ZQj  and  7Q20)  were  determined  by  undertaking  a  single  station 
low  flow  frequency  analysis  for  each  of  the  stations  listed  in  Table  6.1.  Each  prediction 
method  was  then  utilized  in  turn  to  provide  a  simulated  discharge.  Graphical  and  statistical 
(N.S.R.^)  comparisons  were  then  made  for  each  method  as  discussed  in  Section  6.4. 

6.4        Testing  Estimation  Methods 

6.4.1      Testing  of  Regression  Method 

The  preliminary  regression  equations  developed  in  Section  5.2.5  were  tested  for  prediction 
accuracy  using  the  test  stations.  Based  on  these  equations  VQ^  and  7Q20  low  flows  were 
estimated  for  the  test  stations.  The  details  of  testing  are  presented  in  Appendix  A. 

i)  Central  Region 

A  comparison  of  the  observed  and  estimated  low  flows  by  these  regression  equations  is 
shown  in  Figures  6.1(a)  and  6.1(b)  for  7Q2  and  7Q20/  respectively,  along  with  a  45°  line,  on 
which  all  the  pxjints  should  lie  for  f)erfect  simulation.  An  examination  of  these  figures 
indicates  that  except  for  a  few  stations,  the  estimated  flows  are  reasonable  with  an  acceptable 
Nash-Sutdiffe  statistic  (see  Table  6.2). 

ii)  Southeastern  Region 

A  comparison  of  the  observed  and  estimated  low  flows  by  these  regression  equations  is 
shown  in  Figures  6.2(a)  and  6.2(b)  for  7Q2  and  7Q20,  respectively.  These  figures  indicate  an 
imsatisfactory  prediction  for  7Q2  but  satisfactory  prediction  for  7Q20  if  the  two  outlier  stations 
are  excluded  (see  Table  6.2  for  the  Nash-Sutcliffe  statistics). 

The  regression  equations  were  generally  foimd  to  overestimae  both  the  7Q2  and  7Q20  low 
flows  for  most  of  the  stations  where  the  estimation  is  poor. 
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iii)  Discussion  of  Results 

After  a  great  deal  of  effort  and  trying  various  transformations  and  derived  forms  of  various 
parameters,  it  has  been  found  that  the  regression  method  can  lead  to  satisfactory  estimation 
of  low  flows  in  the  Central  region  only.  Ln  the  Southeastern  region,  the  method  gives 
unsatisfactory  estimation  of  7Q2  and  VQ^o  low  flows.  The  development  of  a  regression 
prediction  equation  for  the  combined  regions  did  not  prove  to  be  feasible.  Results  of  the 
testing  are  simtmarized  in  Table  6.2. 

6.4.2      Testing  of  Index  Method 

The  7Q2  low  flows  for  the  Central,  Southeastern  and  the  combined  regions  are  estimated  from 
Figure  5.1  using  the  drainage  area  of  the  respective  station.  The  ZQzo  low  flows  for  all  these 
regions  are  estimated  by  multiplying  the  7Q2  values  by  the  ratio  of  7Q20/7Q2  derived  from 
the  appropriate  figures  for  the  Central,  Southeastern  and  combined  regions.  The  observed 
and  estimated  low  flows  along  with  the  calculation  of  N.S.R.^  for  these  regions  is  presented 
in  App)endix  B.  The  observed  and  simulated  7Q2  and  7Q20  for  the  Central,  Southeastern  and 
combined  regions  are  presented  in  Figures  6.3,  6.4  and  6.5  respectively. 

i)  Central  Region 

Reference  to  Figure  6.3  wiU  indicate  a  reasonable  comparison  of  observed  and  estimated  low 
flows  for  7Q2  and  7Q20.  The  Nash-Sutdiffe  statistic  (see  Table  6.3)  is  also  acceptable. 

ii)  Southeastern  Region 

An  examination  of  Figure  6.4  indicates  that  for  7Q2,  the  flow  is  overpredicted  for  most  of  the 
test  stations.  For  7Q20,  the  flow  is  also  overpredicted  for  a  majority  of  the  stations  but  is 
slightiy  underpredicted  for  02HK003,  Crowe  River  at  Marmora  and  02KC014,  Indian  River 
near  Pembroke.  Station  02HK003  is  regulated  with  a  drainage  area  of  1990  km^  When 
stations  02KC014  and  02LB013  were  excluded  from  the  testing,  the  N.S.R.^  for  7Q20,  improved 
substantially  (see  Table  6.2). 
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iii)  Combined  Regions 

The  results  for  the  combined  regions  indicate  that  the  pattern  of  estimation  for  almost  all  the 
stations  is  the  same  for  7Q20  and  7Q2.  Even  though  the  Southeastern  region  has  fewer  test 
(and  development)  stations  combiriing  the  two  regions  has  not  resulted  in  any  significant 
improvement  in  test  results  of  the  combined  regions  over  the  Southeastern  region  alone. 
When  stations  02KC014  and  02LB013  were  excluded  from  the  analysis,  the  N.S.R.^  for  7Q20 
increased  substantially  (see  Table  6.2). 

iv)   Discussion  of  Results 

It  may  be  generally  concluded  that  the  index  method  is  an  appropriate  method  for  estimating 
7Q2  and  7Q20  low  flows  for  the  Central  region  and  7Q20  in  the  Southeastern  region.  The 
method  is  somewhat  less  reliable  for  both  7Q2  and  7Q20  using  the  result  of  the  combined 
regions.  Although  the  number  of  test  stations  in  the  Southeastern  region  is  fewer  than  the 
Central  region,  combining  the  two  regions  together  did  not  improve  the  prediction  accuracy 
for  the  combined  region  over  the  Southeastern  region. 

6.4.3      Testing  of  Isoline  Method 

The  IsoUne  maps  were  developed  for  the  observed  unit  area  low  flows  of  7Q2  for  the 
development  stations  of  the  Central,  Southeastern  and  the  combined  region  (see  Figures  5.3  a, 
5.3  b,  5.3  c  and  5.3  d).  The  7Q2  unit  low  flows  were  estimated  for  the  test  stations  in  the 
Central  and  Southeastern  regions  by  interpolating  the  values  between  the  contour  Unes.  The 
7Q20  low  flows  were  then  determined  by  multiplying  the  estimated  7Q2  and  7Q20/7Q2  ratio 
determined  from  the  appropriate  figures  for  the  Central,  Southeastern  and  combined  regions 
(i.e.  ratios  developed  in  Index  Method).  The  data  base  for  estimation  of  low  flows  and 
calculation  of  N.S.R.^  is  presented  in  Appjendix  D.  The  results  of  this  analysis  are  presented 
in  Figures  6.6,  6.7  and  6.8  for  the  Central,  Southeastern  and  combined  regions,  resjjectively. 
The  N.S.R.^  for  test  stations  by  the  Isoline  method  is  presented  in  Table  6.2. 
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i)  Central  Region 

The  N.S.R.^  for  IQ^  ai^d  7Q20  is  0.95  and  0.53,  respjectively,  which  is  reflected  in  the  results  on 
Figures  6.6(a),  and  6.6(b).  The  prediction  of  7Q2  is  excellent  and  of  IQ^a  is  fairly  satisfactory. 

ii)  Southeastern  Region 

The  N.S.R.^  for  estimation  of  7Q2  and  7Q20  is  0.66  and  0.72  respectively,  when  stations 
02KC014  and  02LB013  are  included  in  the  analysis.  These  values  are  fairly  satisfactory  for 
the  type  of  system  being  dealt  with  in  this  study.  Figure  6.7(a)  indicates  that  the  prediction 
of  7Qi  low  flow  is  satisfactory  for  most  of  the  stations  and  a  slightiy  overprediction  for  a  few 
stations.  For  the  7Q20  low  flow.  Figure  6.7(b)  indicates  a  satisfactory  estimation  for  most  of 
the  stations  except  slightiy  underestimation  for  station  02HK003,  Crowe  River  at  Marmora. 
This  is  a  regulated  stream  with  1990  km^  drainage  area.  Exclusion  of  stations  02KC014  and 
02LB013  from  the  testing  improved  N.S.R.^  for  estimation  of  7Q2  from  0.66  to  0.89. 

iii)  Combined  Centred  and  Southeastern  Region 

The  N.S.R.^  for  the  combined  region  for  7Q2  and  7<^  is  0.92  and  0.71  respectively,  when 
statior^  02KC014  and  02LB013  were  included  in  the  analysis,  which  indicates  a  satisfactory 
estimation  of  the  low  flows.  Figure  6.8  indicates  a  good  prediction  of  7(^2  and  a  satisfactory 
estimation  of  7Q20.  The  prediction  accuracy  slightly  improved  cifter  exclusion  of  stations 
02KC014  and  02LB013. 

iv)  Discussion  of  Results 

It  is  generally  concluded  that  the  Isoline  method  was  found  to  be  the  most  consistent 
technique  for  satisfactorily  simulating  the  7Q2  and  7Q20  low  flows  for  the  Central, 
Southeastern  and /or  combined  regions. 

6.4.4     Testing  of  Station  Proration  Method 

The  method  of  reciprocal  distance  was  used  in  the  proration  estimation,  i.e.  the  distances 
from  the  test  stations  to  nearby  gauges  were  measured  and  the  weighted  average  of  the  unit 


(MS^vUJ)  SMOTd  a3AU3S90 


(3«^^)  SAAOTd  OAtraSSO 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


^1 


\ 


\ 


\ 


\ 


\ 


\ 


\= 


fi\ 


\ 


\ 


\ 


1 


\ 


I 1 1 1 


E 

■5 

c 
0 

tf) 

<D 

a) 

LU 

C 

05 

n 

c 

in 

n 

CO 

■o 

VI 

<s 

a> 

< 

^ 

i> 

<■> 

3 

0 

01 

u. 

5 

o 

0 

g 
"cB  O 

5  E 

I   ^ 

U_     03 

If 
o 


a> 

a> 

c 

i 

'J 

a. 

-0 

ti 

s 

9 

iO 

^ 

^ 

IH 

S 

'5> 

u 

M 

a 

E 

^ 

a 

3 

e 
0 

u 

U 

(oes^^)  SMOTd  (HAH3SQO 


T ^ 3 '  T 

(3««,uj)  SAAOTd  a3AH3S90 


6-6 

low  flows  was  estimated  for  the  test  stations  by  using  the  reciprocal  of  distances  to  nearby 
stations  and  their  observed  unit  area  low  flows.  The  figures  from  which  these  distances  were 
measured  and  the  tables  presenting  the  data  base  for  estimation  of  low  flows  and  calculation 
of  Nash-Sutdiffe  R^  are  presented  in  Appendix  C.  The  N.S.R.^  for  estimation  of  VQ^  and  7Q20 
for  the  Central,  Southeastern  and  combined  regions  are  presented  in  Table  6.3.  The  results 
of  this  analysis  for  these  regions  are  presented  in  Figures  6.9,  6.10  and  6.11. 

i)  Central  Region 

The  prediction  of  7Q2  low  flow  is  satisfactory  as  is  evidenced  by  Figure  6.9(a)  while  for  7Q20 
the  flow  for  02EB004,  North  Branch  Muskoka  River  at  Port  Sydney  is  overpredicted  although 
within  acceptable  limits  for  all  other  stations.  (Station  02EB004  is  a  regulated  stream  with  a 
drainage  area  of  1390  km^). 

ii)  Southeastern  Region 

An  analysis  of  Figures  6.10(a)  and  6.10(b)  indicates  that  both  the  7Q2  and  7Q20  low  flows  are 
overpredicted  for  almost  all  the  stations  except  02HK003,  Crowe  River  at  Marmora.  The 
prediction  for  this  region  by  this  method  is  unsatisfactory  for  both  the  7Q2  and  7Q20  low 
flows.  The  exclusion  of  stations  02KC014  and  02LK013  did  not  improve  the  result 
appreciably. 

iii)  Combined  Central  and  Southeastern  Regions 

The  N5.R.^  for  7Q2  for  the  combined  region  is  acceptable  but  not  for  7Q20.  An  examination 
of  Figures  6.11(a)  and  6.11(b)  indicates  that  the  prediction  for  7Q2  is  satisfactory.  However, 
for  7Q20,  the  results  are  not  as  good  no  matter  if  the  analysis  is  conducted  by  including  or 
excluding  the  stations  02KC014  and  02LK013. 


6.5        Summary 


The  Nash-Sutcliffe,  R^  for  the  test  stations  for  all  the  regions  and  different  methods  of 
estimation  of  low  flows  is  given  in  Table  6.2.  The  results  are  summarized  below  for  each 
region. 
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i)  Central  Region 

The  N.S.R.^  for  estimation  of  70:  is  the  highest  (0.96)  by  the  station  proration  method 
followed  by  the  isoline  method  (0.95),  index  method  (0.83)  and  the  regression  method  0.81). 
For  7(^20  the  highest  (0.84)  is  by  the  regression  method  followed  by  the  index  method  (0.82), 
isoline  method  (0.53)  and  station  proration  (0.02).  It  is  concluded  that  the  index,  regression 
or  isoline  methods  may  be  used  to  estimate  low  flows  with  acceptable  accuracy  for  this 
region. 

ii)  Southeastern  Region 

The  N.S.R.^  for  estimation  of  7Q2  is  the  highest  (0.89)  for  the  Isoline  method  followed  by 
regression  (0.87),  station  proration  (0.12)  and  Index  method  (-0.12),  while  for  7(^20  the  highest 
(0.92)  is  by  the  Index  method  followed  by  Isoline  method  (0.73),  regression  (0.84)  and  station 
proration  (-0.18).  It  is  concluded  that  the  Isoline  method  can  satisfactorily  simulate  both  7Q2 
and  7Q20  in  this  region.  Since  the  7Q20  low  flows  are  estimated  from  the  maps  of  7Q2  and 
multiplying  it  by  a  VQ^o/VQj  ratio  as  described  in  the  Index  method,  it  can  be  anticipated  that 
the  other  low  flows  can  also  be  satisfactorily  simulated  by  these  methods. 

iii)  Combined  Central  and  Southeastern  Regions 

The  N.S.R.^  for  estimation  of  7Q2  is  highest  (0.95)  by  the  Isoline  method  followed  by  station 
proration  (0.86)  and  Index  method  (0.45).  For  7Q20  the  best  result  (0.75)  is  by  the  Isoline 
method,  followed  by  the  index  method  (0.44)  and  station  proration  (0.04).  A  satisfactory 
regression  equation  could  not  be  produced  when  the  regions  were  combined.  Oi\Iy  the 
isoline  method  could  satisfactorily  estimate  the  low  flows  for  the  combined  Central  and 
Southeastern  regioits. 

iv)  Discussion  of  Results 

For  estimation  of  low  flows  in  the  Central  region,  any  one  of  the  index,  regression  or  isoline 
could  be  used  satisfactorily.  On  the  other  hand,  for  the  Southeastern  and  the  combined 
regions,  only  the  isoline  method  produced  cor^istent  estimates  of  low  flows.  It  is,  therefore. 
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concluded  that  the  isoline  method  is  the  most  robust  technique  for  estimating  low  flows  in 
the  Central  and  Southeastern  regions. 

It  should  be  noted  that  all  the  estimation  technqiues  could  not  produce  satisfactory  results 
for  some  stations.  Therefore,  care  should  be  taken  when  using  these  methodologies  in 
estimating  low  flows  for  ungauged  watersheds. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 

7.1  Conclusions 

i)   Central  Region 

a)  For  this  region,  any  of  the  methods,  i.e.  Index,  Regression  or  Mapped  Isoline  can  be 
used  for  estimating  low  flows. 

b)  The  regression  equation  developed  is  based  on  the  parameters  DA  and  BFI. 

c)  The  regression  equations  based  on  these  parameters  were  developed  for  the  3,  7  and 
30  day  low  flows  for  recurrence  intervals  of  2, 20  and  50  years  and  presented  in  Table 
5.4. 

ii)  Southeastern  Region 

a)  Three  of  the  estimation  techniques  (Regression,  Index  and  Station  Proration)  did  not 
provide  coi^sistently  robust  estimates  of  the  7Q2  and  7Q20  low  flows.  The  Index 
method  gave  excellent  prediction  for  7Q20  but  very  poor  for  7Q2.  For  the  regression 
method,  the  equation  based  on  DA^  and  BFI  gave  a  satisfactory  prediction  for  7Q2  but 
unsatisfactory  for  VQ^q.  It  is  concluded  that  none  of  these  methods  (Index,  Regression 
or  Station  Proration)  could  satisfactorily  simulate  the  low  flows  for  this  region. 

b)  The  Mapped  Isoline  method  provided  the  most  consistently  robust  method  for 
estimation  of  both  the  7Q2  and  7Q20  low  flows. 

iii)  Combined  Central  and  Southeastern  Regions 

a)  The  Index  Method  resulted  in  fair  estimation  of  low  flows  while  the  station  proration 
method  is  not  consistent  in  estimating  both  the  7Q2  and  7Q20  low  flows,  resulting  in 
a  satisfactory  prediction  for  7Q2  and  poor  for  7Q20.  A  satisfactory  regression  equation 
could  not  be  produced  for  the  combined  regions. 
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b)       The  Mapped  Isolines  method  gave  good  prediction  for  7Q^  and  7Q2o- 

Overall,  it  was  concluded  that  the  Mapped  Isoline  Method  was  the  most  consistent  for 
prediction  of  both  the  7Q2  and  7Q20  low  flows. 

7.2         Recommendations 

a)  Isoline  maps  of  unit  area  low  flows  should  be  prepared  by  digital  terrain  modelling 
using  a  Geographical  Information  System  (GIS)  for  all  the  regions  defined  by  the 
Ministry  of  the  Environment. 

b)  Refinement  of  the  regression  technique  for  watersheds  in  the  Southeastern  region 
appears  to  require  additional  investigation  including  definition  and  incorporation  of 
other  watershed  parameters  in  order  to  suffidentiy  explain  low  flow  characteristics. 

c)  The  regression  studies  suggest  that  the  maximum  fluctuation  of  the  groundwater 
table  may  be  an  important  parameter  under  certain  conditions.  The  maximxmi 
fluctuations  of  the  water  table  for  a  well  was  calculated  from  the  observed  maximum 
and  minimum  water  table  depths  repwrted  in  the  Environment  Adas  of  Observation 
Wells  in  Ontario  published  by  the  Ontario  Ministry  of  the  Environment,  1980.  The 
data  was  old  and  only  sparse  well  observation  information  was  available,  so  it  is 
recommended  that  the  Atlas  of  the  Observation  Wells  be  updated  and  records  for  a 
larger  number  of  observation  wells  should  be  obtained. 

d)  Statistical  tests  indicate  that  a  large  number  of  the  stations  data  used  in  this  study 
show;  trend,  dependence  and  some  non-randomness  characteristics.  Future 
investigations  should  examine  and  quantify  the  extent  of  trends  on  low  flow  series. 

e)  The  single  station  analysis  was  based  on  flow  data  up  to  1986.  Since  5  years  of 
additional  data  for  active  stations,  and  because  some  additional  stations  record  length 
may  now  qualify  for  reliable  statistical  analysis,  we  feel  the  single  station  data  should 
be  updated  for  future  studies  and  regular  analysis  of  gauge  station  data. 
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The  techniques  developed  within  this  study  do  not  always  accurately  depict  low 
flows  for  all  stations  and,  therefore,  care  is  recommended  by  anyone  applying  these 
techiuques. 
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APPENDIX  A 
REGRESSION 

A.1        Development  of  Regression  Equations 

Over  250  preliminary  regression  equations  were  developed  in  order  to  predict  the  7Q2  and 
7Q20  low  flows  as  a  function  of  hydrometeorologic  and  basin  physiographic  parameters, 
parameter  transforms  or  derived  transforms.  Potential  equations  were  evaluated  based  on 
various  groups  of  parameters  (see  Table  A.l). 

Selected  preliminary  regression  equations  were  developed  and  tested  for  estimation  of  the 
7Q^  and  7Q20  low  flows  in  the  Central,  Southeastern  and  combined  regioris. 

On  the  basis  of  the  highest  value  of  Nash-Sutcliffe  R^,  for  the  testing  stations,  to  be  described 
subsequently,  the  following  groups  of  parameters  were  found  to  give  the  highest  prediction 
accuracy: 

Central  Region:  DA  and  LOGjo  EVA 

Southeastern  Region:  DA^  BFP  and  LNTH^ 

Combined  Regions:  DA^  BFI^  and  LNTH^ 

The  coefficients  of  these  parameters  for  different  low  flows  are  presented  in  Tables  A.l,  A.2 
and  A.3  for  the  Central,  Southeastern  and  the  combined  region.  These  results  are  obtained 
when  two  stations,  02KC014  (Iitdian  River  near  Pembroke)  and  02LB013  (South  Nation  River 
near  Cassehnan)  were  included  in  the  test  stations  for  the  Southeastern  region  and  hence  in 
the  combined  region.  When  these  stations  were  excluded,  then  for  the  combined  region,  the 
regression  equation  based  on  DA^,  BFI^  and  Logjo  EVA  gave  the  best  estimation  of  both  the 
7Q2  and  7Q20.  The  coefficients  of  the  regression  equations  based  on  these  parameters  for  the 
combined  region  are  given  in  Table  A4. 


Table  A.  1  Summary  of  Regression  Analysis 


GrouD 

Repon 

Parameter/Transformed  or  Denved  Forrr. 

Low  Flow 

N 

Group  I 

CN 

MAP,  DA.  LNTH.  RN.  GW,  BR.  MAS.  MAE.  SLP.  ACLS  and  MAR 
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•18 

7O20 

Jg 

SE 

702 

25 

7O20 

25 

CNSE 

702 

-•< 

t7O20 

■3 

jrou3  II 

1 

Fint 
Set 

CN 

CW^3.  R_N.  MAP,  LoelO  MA£.  SLP'3.  SF2  ANT)  LoelO  BFI 

"702 

48 

LoElO  GW,  R,N,  INT,  MAS,  SLP':,  STZ.  BFl'3  and  ACLS 

7O20 

48 

SE 

LoBlO  GW,  R.N,  INF,  SF3,  BFr3  and  ,^1115 

■702 

25 

LoelO  GW.  R,N,  INT,  SF3,  BF!'3  and  .ACLS 

■7O20 

2*; 

CNSE 

GW^3.  R.N.  M.\P'-i.  MAR''3,  MAE,  ,MAS^3,  SF?,  BR"}  and  ACLS 

702 

--, 

GWM,  R>-.  .MAP"3,  LoRlO  .MAR,  .MAE'3,  .VUS"3,  SF3  ANT)  BR^i 

7O20 

"3 

Second 
Set 

CN- 
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73 

Group  III 

CN 

MAE  and  DA 
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48 
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SE 
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25 
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25 

CNSE 

la)  DA.  GW  and  MAE 

[702 

73 

7O20 

73 

(b)  DA.  BFI  aod  MAE 

702 

73 

7O20 

73 

Group  IV 

CN-^-^^-t    - 

OA  aDiI,^mjSiU£^^,i-,^-.iiiitfi— _i:^T».rc-  »  -, rr-sa=s-s9r^-. 

701. — 

-      4« 

7O20--- 

■-  4« 

SE 

DA'3.  MAE'3  and  Bn'2 

702 

25 
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25 

CNSE-  _. ,  - 

(a)  DA'3.  Log  10  MAE  and  LoglO  GW 

702 

73 

7O20 

73 

[b)  DA«3,  LoglO  MAE  and  BH-^S 

702 

73 

7020 

73 

Group  V 

CN 

DA-^V  BIT^iand  IJ<TH*2 

702 
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7O20 

48 

SE 

702--- 

--    25 

7O20- 

25 

CNSE 

702 

73 

7O20 

73 

Group  VI 

CN 

DA  only 

702 

48 

7020 

48 

SE 

702 

25 

7O20 

25 

CNSE 

702 

73 

7O20 

73 

N-Number  of  Stations  SEE-Standard  Error  of  Estimates 


Table  A-2  Summary  of  Regression  Analysis  (Central  Region) 
Y-aO+al*DA+a2*LOOIO  EVA 


Dependent 
Parameter 

Independent  Parameters 

aO 

al 

a2 

N 

SE 

R'>2 

7O20 

27.32 

3.02*10^-4 

-9.33 

48 

0.29 

0.62 

702 

130.24 

2.39*  10^-4 

-44.69 

48 

0.624 

0.8 

302 

50.42 

9.58*10^-4 

-17.27 

48 

0552 

072 

3O20 

22.51 

2.94*10^-4 

-7.68 

48 

0.266 

06 

3O50 

17.36 

2.39*10^-4 

-5.91 

48 

0.241 

0.53 

30O2 

75.7 

1.24*10^-3 

-25.95 

48 

0.764 

0.72 

30O20 

39.23 

4.39*10^-4 

-13.42 

48 

0.408 

0.63 

30O50 

33.48 

4.03*10^-4 

-11.45 

48 

0374 

061 

Table  A-  3 Summary  of  Regression  Analysis  (South  Eastern  Region) 
Y=aO+a  1  •DA^3+a2*Bn''2+a3*LNTH^2 


Dependent 
Parameter 

Independent  Parameters 

aO 

al 

a2 

a3 

N 

SE 

R^2 

7O20 

-0529 

7.31*10^11 

1.759 

2.31*10^-5 

48 

091 

087 

702 

-0.78 

1.24*10^-10 

3.186 

3.53*10^5 

48 

1.37 

089 

302 

-0.95 

1.24*10^-10 

3.935 

-2.45*10^-5 

48 

1.25 

091 

3O20 

-0.73 

7.37*10^-11 

2.668 

-3.87*10^-5 

48 

0.84 

088 

3O50 

-0.78 

6.86*10^-11 

2.877 

-5.01*10^-5 

48 

081 

088 

30O2 

-0.76 

1.39*10^-10 

3.289 

4.74*10^-5 

48 

0.9 

09 

30O20 

-0.38 

7.79*10^-11 

1.102 

7.68*10^-5 

48 

1.1 

0.84 

30O50 

-033 

7.18*10^11 

0869  18.31*10^-5 

48 

1.047 

083 

Table  A -4  Summary  of  Regression  Analysis  (Combined  Central  and  South  Eastern  Region) 
Y=aO+a  1  *DA'^3+a2*B^'^2-^a3*LNTH'^2 


Dependent 
Parameter 

Independent  Parameters 

aO 

al 

a2 

a3 

N 

SE 

R^2 

7O20 

-0215 

7.32*10^-11 

1.887 

-1.55*10^-5 

73 

0.649 

087 

702 

-0666 

1.16*10^10 

5.254 

-5.80*10^-5 

73 

1.23 

085 

302 

-0.641 

1.17*10^10 

4.393 

-4.11*10^-5 

73 

098 

09 

3020 

-0.287 

7.25*10^-11 

2.03 

-3.46*10^-5 

73 

061 

088 

3O50 

-0.278 

6.85*10^-11 

1.785 

-3.00*10^-5 

73 

058 

088 

30O2 

-0.688 

1.28*10^-10 

5.35 

-2.33*10^-5 

73 

1.24 

088 

30O20 

-0.239 

7.70*10^11 

2.38 

-1.04*10^-5 

73 

0.81 

0.84 

30O50 

-0.206 

7.19*10^-11 

2.05 

-3.29*10^-5 

73 

0.75 

084 

Table  A-5    Summary  of  Regression  Analysis  (Combined  Central  and  South  Eastern  Region) 
Y=a0+al*DA'^3■^a2*B^^3-(■a3*LOG10  EVA 


Dependent 
Parameter 

Independent  PaiM^eters 

aO 

al 

a2 

a3 

N 

SE 

R'>2 

7O20 

40.08 

6.35*10^-11 

1.00 

-13.82                     73 

0.57 

0.90 

702 

105.36 

8.94*10^-11 

2.69 

-36.35 

73 

0.92 

092 

302 

62.29 

9.73*10^-11 

3.23 

-21.55 

73 

083 

093 

3O20 

33.15 

6.13*10^-11 

1.26 

-11.47 

73 

055 

0.90 

3O50 

28.46 

5.86*10^-11 

1.17 

-9.86 

73 

053 

0.90 

30O2 

86.13 

1.06*10^-10 

3.79 

-29.73              73 

1.02 

092 

30O20 

53.64 

6.55*10^-11 

1.20 

-18.48              73 

0  69 

089 

30O50 

49.20 

6.21*10^-11 

1.00 

-16.95              73 

0.65 

088 

A-2 
A^       Testing  of  the  Regression  Equations 

On  the  basis  of  the  testing  and  development  of  the  regression  equations,  engineering 
judgement  and  the  practical  utility  of  these  equations,  the  following  groups  of  parameters 
were  found  to  have  given  the  most  satisfactory  results: 

Central  Region:  DA  and  BFI 

Southeastern  Region:  DA^  and  BFI 

It  was  found  that  no  common  regression  equation  could  satisfactorily  predict  the  low  flows 
for  both  the  regions.  This  indicated  that  the  prediction  of  the  low  flows  by  regression  should 
be  carried  out  independently  for  each  region. 

The  coefficients  of  those  parameters  for  different  low  flows  are  presented  in  Tables  5.4  and 
5.5  in  the  main  text,  for  the  Central  and  Southeastern  regions,  respectively.  Since  the 
prediction  of  low  flows  by  regression  method  is  not  consistent  with  the  development  of 
regression  equation  for  the  Southeastern  region  as  regards  the  accuracy  of  prediction  is 
concerned,  it  was  tried  to  simulate  7Q20  from  7Q2  and  using  the  ratio  of  7Q20/7Q2  from  Figure 
5.2  as  0.228.  The  7Q2  low  flows  were  predicted  by  the  regression  equation.  This  method 
improved  the  accuracy  of  prediction  of  7Q20  tremendously  over  the  prediction  of  7Q20  directly 
by  regression  equation. 

Further  investigations  indicated  that  when  two  heavily  regulated  stations  (02KC014  and 
02LB013)  were  removed  from  the  testing,  the  regression  method  produced  satisfactory  results 
for  the  7Q20  estimates  for  the  test  stations  in  the  Southeastern  region. 


TABLE  A.6 
SUMMARY  OF  HNAL  REGRESSION  ANALYSIS  TESTING 


Group 

Parameters 

NASH-SUTCLIFFE  R' 

Region 

Central 

Southeastern 

III 

EVA,  DA,  BR  and  GW 

0.71 

0.72 

-6.32 

-15.49 

IV 

DA^  BFP  and  LNTH" 

028 

0.4 

-0.66 

-0.51 

V 

DA 

0.68 

0.67 

-7.35 

-11.32 

VI 

SF2,  MAP,  MAR,  EVA  and  BFI 

0.82 

0.78 

-443 

-6.81 

vn 

SF2  and  BFI 

0.65 

0.72 

-4.3 

-6.8 

DC 

DA,  LNTH  and  BFI 

0.8 

0.85 

-7.3 

-11.2 

X 

DA  and  BFI 

0.81 

0.84 

-7.1 

-11.9 

XI 

DA^  and  BH 

0.76 

0.75 

-1.4 

-3.1 

xn 

DA"  and  BFI 
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0.84 
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-16.1 
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APPENDIX  B 
INDEX  METHOD 

The  observed  and  estimated  low  flows  by  Index  Method  along  with  the  calculation 
of  Nash-Sutdiife,  R^,  for  the  Central,  Southeastern  and  combined  regions  are 
presented  in  Tables  B.l,  B.2  and  B.3,  respectively. 

The  7Q2  low  flows  for  the  Central,  Southeastern  and  the  combined  regions  are 
estimated  from  Figure  5.1  using  the  drainage  area  of  the  respective  station.  The  7Q20 
low  flows  for  all  these  regions  are  estimated  by  multiplying  the  ZQj  values  by  the 
ratio  of  7Q20/7Q2  derived  from  the  appropriate  figures  for  the  Central,  Southeastern 
and  combined  regions.  The  observed  and  estimated  low  flows  along  with  the 
calculation  of  N.S.R.^  for  these  regions  is  presented  in  Appendix  B.  The  observed  and 
simulated  7Q2  and  7Q20  for  the  Central,  Southeastern  and  combined  regions  are 
presented  in  Figures  6.3,  6.4  and  6.5  respectively. 
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APPENDIX  C 
ISOLINE  METHOD 

The  data  base  for  estiinating  low  flows  by  isoline  method  is  presented  in  Tables 
D.l(a)  and  D.2(a)  and  D.3(a)  for  the  Central,  Southeastern  and  combined  regions, 
respectively.  The  summaries  of  the  simulation  results  by  this  method  along  with  the 
calculation  of  Nash-Sutdiffe,  R'',  for  the  Central,  Southeastern  and  the  combined 
regions  are  presented  in  Tables  C.l,  C.2,  and  C.3  respectively. 

The  Isoline  maps  were  developed  for  the  observed  unit  area  low  flows  of  JQ^  for  the 
development  stations  of  the  Central,  Southeastern  and  the  combined  region  (see 
Figures  5.4, 5.5  and  5.6).  The  70^  unit  low  flows  were  estimated  for  the  test  stations 
in  the  Central  and  Southeastern  regions  by  interpolating  the  values  between  the 
contour  lines.  The  VQ^q  low  flows  were  then  determined  by  multiplying  the  estimated 
7Q2  and  7Q20/7Q2  ratio  determined  from  the  appropriate  figures  for  the  Central, 
Southeastern  and  combined  regions  (i.e.  ratios  developed  in  Index  Method).  The  data 
base  for  estimation  of  low  flows  and  calculation  of  N.S.R.^  is  presented  in  Appendix 
D.  The  results  of  this  analysis  are  presented  in  Figures  6.6,  6.7  and  6.8  for  the 
Central,  Southeastern  and  combined  regions,  respectively.  The  N.S.R.^  for  test  stations 
by  the  Isoline  method  is  presented  in  Table  6.2. 
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APPENDIX  D 
STATION  PRORATION  METHOD 

The  data  base  for  estiinating  low  flows  by  the  station  proration  method  is  presented 
in  Tables  D.l  cind  D.2  for  the  Central  and  Southeastern  regions,  respectively.  The 
summaries  of  the  simulation  results  by  this  method  along  with  the  calculation  of 
Nash-Sutdiffe,  R^  for  the  Central,  Southeastern  and  the  combined  regions  are 
presented  in  Tables  D.3,  D.4  and  D.5  respectively. 

The  method  of  reciprocal  distance  was  used  in  the  proration  estimation,  i.e.  the 
distances  from  the  test  stations  to  nearby  gauges  were  measured  and  the  weighted 
average  of  the  imit  low  flows  was  estimated  for  the  test  stations  by  using  the 
reciprocal  of  distances  to  nearby  stations  and  their  observed  unit  area  low  flows.  The 
figures  from  which  these  distances  were  measured  and  the  tables  presenting  the  data 
base  for  estimation  of  low  flows  and  calculation  of  Nash-Sutdiffe  R^  are  presented  in 
Appendix  C.  The  N.S.R.^  for  estimation  of  7Q2  and  ZQzq  for  the  Central,  Southeastern 
cmd  combined  regions  are  presented  in  Table  6.3.  The  results  of  this  analysis  for  these 
regions  are  presented  in  Figures  6.9,  6.10  and  6.11. 
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E-1 

APPENDIX  E 

ANALYSES  OF  EXISTING  TECHNIQUES 


To  assess  the  applicability  of  the  previous  results  extrapolated  to  the  present  study 
area,  the  previously  develof)ed  techniques  wdth  the  previously  determined  coefficients 
and  constants  were  utilized  to  predict  low  flow  characteristics  for  statior\s  located  in 
the  Central,  Southeastern  Regions  and  the  combined  regions.  The  analysis  is 
summarized  in  Tables  E.l,  E.2  and  E.3  for  the  above  regions  respectively. 
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TABLE  E^ 

TESTING  PREVIOUS  REGIONAL  ANALYSIS 

OF  LOW  FLOW  CHARACTTERISTICS 

IN  THE  SOUTHEASTERN  REGION 
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F-1 
APPENDIX  F 

SKVIPLE  CORRELATION  OF  PARAMETERS 

The   following   Table   F.l    is   a   matrix   of   the   correlation   coefficients   for   the 
untransformed  independent  parameters  used  in  this  study. 


TABLE  F.1 

SIMPLE  CORRELATION  OF  PARAMETERS 

CE^JTRAL  REGION 


CORRELATIONS  /VSRIOBLES  SK  TO  020. 


nAP  nes 

.2186  -.0308 

-.0220  -Am 
l.OODO     .6877t: 

.687711  1.0000 

.U19H  .1567 

-.37171  -.4«8«1I 

-.S124JI  -.626811 
.510011  .M07I1 
.«S2II  .:4('BI 
.34011  .520011 
.36651  .33681 
.610611  .461111 
.474611    .469611 


H  ot  cases:   46       l-tailei  Sionit:  I  -  .01  11  -  .001 
'  .  ■  i;  cnnteJ  if  a  :5e^'i:ient  cannot  be  cojputei) 


eiations:  iU 

RN 

1.0000 

.1252 

.1252 

1.0000 

"       .2186 

-.0220 

5       -.0308 

-.1908 

R       .0268 

.1135 

F       -.096? 

-.2146 

a       -.1395 

.1062 

.0466 

-.1160 

!        .0314 

.)1545 

TH      -.0391 

-.1264 

LS       .0592 

.1633 

.0000 

.0512 

.555011 

.2805 

.437411 

.2248 

.0800 

.463711 

.42741 


SLP 

-.0969 
-.2146 
-.37171 
-.448411 

.0512 
1.0000 

.3273 
-.38311 
-.0691 
-.39161 
-.42571 
-.1544 
-.0650 


EVA 

-.1395 

.1062 
-.892611 
-.626811 
-.555011 

.3273  " 
KOOOO 
-.548811 
-.485011 
-.36641 
-.41861 
-.710611 
-.497411 


DA 

.0466 

-.1160 
.510011 
.540711 
.2305 

-.38311 

-.548811 

1.0000 
.33951 
.813011 
.36061 
.714711 
.662111 


.0814 

.0545 

.486211 

.34061 

.637411 
-.0691 
-.485011 

.33951 
1.0000 

.2485 

.2207 

.446511 

.456911 


LNIH 

-.039! 
-.1264 

.34011 

.520011 

.2248 
-.39181 
-.36641 

.813011 

.2485 
1.0000 

.2610 

.445311 

.474311 


Correlations:  ACLS 

EV  .0592 

RN  .1633 

BSP  .36651 

BAS  .33681 

flOR  .0800 

SLR  -.42571 

EVA  -.41861 

DA  .36061 


LNIH 
ACLS 
02 


.2610 
1.0000 
.2497 


.0850 
-.1835 
.610611 
.161111 

.4o3:'ii 

-.1544 

-.710611 
.714711 
.446511 
.4(5311 

.2697 
1.0000 
.882511 


.0277 

.47)«ll 

.469611 

.42741 

.0650 

.497411 

.662111 

.456911 

.474311 

!3:;5ii 


SOUTHEASTERN  REGION 


C:SREL«TIONS  .'VARIABLES  GW  TO  02i 


Correlations: 

ED 

RN 

611 

1.0000 

.2097 

M 

.2097 

1.0000 

•5P 

.56151 

.2034 

riAS 

-.54421 

-.4268 

m 

.3537 

-.1846 

SLP 

-.0007 

-.0598 

EVA 

.654611 

.2629 

DA 

-.1956 

.3512 

EH 

.278!. 

.744811 

LIITH 

.0127 

.4444 

S-L5 

.4295 

.773411 

02 

-.2543 

.2999 

020 

-.2756 

.2732 

HAP 

HAS 

HAR 

SLP 

EVA 

DA 

BFI 

LNIH 

Correlations:  ACLS 

02 

020 

..56151 

-.54421 

.3537 

-.0007 

.654611 

-.1956 

.2783 

.0127 

EM 

.4295 

-.2543 

-.2754 

.2034 

-.4268 

-.1866 

-.0598 

.2629 

.3512 

.744811 

.4444 

RN 

.773411 

.2999 

.2782 

1.0000 

-.639911 

.4068 

.0913 

.47851 

-.3045 

.2731 

-.2364 

nur 

.3149 

-.2837 

-.2992 

-.639911 

1.0000 

-.1315 

-.1023 

-.50931 

.0352 

-.3977 

-.0779 

riAS 

-.50011 

.1174 

.1-56 

.4048 

-.1315 

l.OOOO 

-.55441 

.52771 

-.3392 

-.1879 

-.2272 

BAR 

-.0578 

-.2710 

-.2464 

.0913 

-.1028 

-.55441 

1 .  0000 

-.48101 

-.1345 

.1219 

-.3049 

SLP 

.1475 

-.0573 

-.0491 

.47851 

-.50931 

.52771 

-.48101 

1.0000 

-.2760 

-.0026 

-.0364 

EVA 

.2452 

-.4544 

-.3045 

.0352 

-.3392 

-.1345 

-.2740 

1.0000 

.4455 

.773911 

DO 

.0739 

.900711 

.887711 

.2731 

-.3977 

-.1879 

.1219 

-.0026 

.4455 

1.0000 

.51221 

PFl 

.642911 

.4546 

.4325 

-.2364 

-.0779 

-.2272 

-.3049 

-.0366 

.773911 

.51221 

1.0000 

LNTH 

.2047 

.430211 

.3149 

-.50011 

-.0578 

.1475 

.2452 

.0739 

.642911 

.2047 

ACLS 

1.0000 

.0990 

.0781 

-.2837 

.1174 

-.2710 

-.0573 

-.4564 

.900711 

.4568 

.430211 

02 

.0990 

l.OnOO 

.997411 

-.2992 

.1388 

-.2664 

-.0491 

-.46421 

.887711 

.4325 

.e21Bll 

020 

.0781 

.997411 

l.OOOO 

1  of  cases:       25  l-tailed  Signif:    I  -  .01    11 

.  '  ;=  printed  if  a  coefficient  cannot  be  coiputed 


COMBINED  REGIONS 


C2RFELATI0NS  /vaPIt3LES  C-»  TO  020. 


SPSS/PC* 

Correlations:  3W 

RN 

BAP 

BAS 

IMR 

SLP      EVA 

DA 

BFI 

LNIH 

Correlations:  ACLS 

2 

02" 

5.1 

l.OOOO 

.1537 

.23591 

-.0424 

.2430 

-.1772 

0324 

.0288 

.1768 

.1740 

E« 

.375711   - 

0957 

-.1331 

^^' 

.1537 

1.0000 

.0334 

-.1893 

-.0224 

-.1601 

1112 

.1576 

.31211 

.1448 

RN 

.370311 

0377 

.0473 

r.if 

.2S591 

.0334 

1.0000 

.579011 

.432211 

-.2413 

462211 

.0447 

.413511 

.1197 

BAP 

.33701 

1616 

.0402 

HAS 

-.0424 

-.1893 

.577011 

l.OOOO 

.1505 

-.359611   - 

•92411 

.2133 

.2181 

.30591 

BAS 

.2258 

2371 

.1392 

f.'riy 

.2430 

-.0224 

.432211 

.1505 

1.0000 

-.1845 

370911 

.0430 

.389011 

.1570 

BAR 

.1352 

1424 

.0685 

SLP 

-.1772 

-.1401 

-.2413 

-.359411 

-.1845 

1.0000 

1917 

-.27651 

.0010 

-.422911 

SLP 

-.29021 

1271 

EVA 

-.0324 

.1142 

-.442211 

-.492411 

-.370911 

.1917    I 

0000 

-.436611 

-.2571 

-.374511 

EVA 

-.32561   - 

-.413011 

DA 

.0238 

.1573 

.0o47 

.2133 

.0630 

-.27451   - 

436411 

1.0000 

.34631 

.785611 

DA 

.2431 

6507II 

.815311 

F-I 

.1766 

.31211 

.413511 

.2131 

.339011 

.0010 

2571 

.34631 

1.0000 

.35351 

BFI 

.372511 

418011 

.361111 

LNTH 

.1740 

.1448 

.1197 

.30591 

.1570 

-.422911   - 

374511 

.785611 

.35351 

l.nOOO 

LNIH 

.33601 

572111 

;CLE 

.375711 

.370311 

.33701 

.2258 

.1352 

-.29321 

32541 

.2431 

.372511 

.33401 

ACLS 

1.0000 

1904 

.1213 

32 

-.0957 

.0877 

.1416 

.2371 

.1426 

-.1271 

495311 

.350711 

.416011 

.572111 

02 

.1906    1 

0000 

020 

-.1331 

.0473 

.0402 

.1892 

.0685 

-.0859 

413011 

.845311 

.381111 

.564311 

820 

.1213 

970911 

1.0000 

H  cf  cases:        "  l-tailed  Signif:     I  -  .01     11 

■   .   ■  IS  jrnted  :f  a  coefficient  canrot  be  cciouted 
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